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CONVERSION  FACTORS,  NON-SI  TO  SI  (METRIC) 
UNITS  OF  MEASUREMENT 


Non-SI  units  of  measurement  used  in  this  report  can  be  converted  to  SI  (metric) 
units  as  follows: 


Multiply 

By 

To  Obtain 

degrees  (angle) 

0.01745329 

radians 

feet 

0.3048 

metres 

horsepower  (550  foot-pounds 
(force)  per  second  per 
ton  (force) 

83.82 

watts  per  kilonewton 

inches 

2.54 

centimetres 

miles  (US  statute) 

1.609347 

kilometres 

pounds  (force) 

4.448222 

newtons 

pounds  (force) -inch-second 
squared 

0.11306064 

kilograms -square  metres 

pounds  (force)  per  square  inch 

6.894757 

kilopascals 

pounds  (mass) 

0.4535924 

kilograms 

pounds  (mass)  per  cubic  foot 

16.01846 

kilograms  per  cubic  metre 

pounds  (mass)  per  cubic  inch 

0.0276799 

kilograms  per  cubic  centimetres 

square  inches 

6.4516 

square  centimetres 

STEERABILITY  ANALYSIS  OF  TRACKED  VEHICLES; 
THEORY  AND  USER'S  GUIDE  FOR  COMPUTER  PROGRAM  TVSTEER 


CHAPTER  1 
INTRODUCTION 


1.1  BACKGROUND 

A  terrain-tracked  vehicle  interaction  model  for  predicting  the  steering 
performance  of  ground-crawling  vehicles,  referred  to  as  TVSTEER,  was  developed 
in  1978  (Reference  1)  and  validated  in  1981  and  1982  (References  2  and  3). 

The  model  treats  both  uniform  (steady-state)  and  nonuniform  (transient) 
motion. 

In  1985,  the  Technical  Management  Committee  (TMC )  for  the  NATO  Reference 
Mobility  Model  (NRMM)  recommended  the  acceptance  of  TVSTEER  as  a  mobility 
support  model  for  evaluating  high-mobility/agility  tracked  combat  vehicles. 

The  TMC  also  recommended  that  a  TVSTEER  user's  guide  be  prepared  by  the  US 
Army  Engineer  Waterways  Experiment  Station  (WES)  for  the  NATO  community. 

1.2  PURPOSE 

This  report  summarizes  the  mathematical  development  of  the  terrain- 
tracked  vehicle  interaction  model  and  presents  a  user's  guide  for  the  computer 
program  TVSTEER. 

1.3  SCOPE 

The  development  of  the  terrain-vehicle  interaction  model  is  given  in 
Chapter  2.  A  user's  guide  for  the  computer  program  TVSTEER,  including  typical 
input  and  output  data,  flowcharts  of  the  main  program  and  two  subroutines,  and 
several  sample  runs,  is  presented  in  Chapter  3-  Appendix  A  presents  the 
program  listing  for  TVSTEER,  and  Appendix  B  is  a  glossary  of  important 
variables  for  the  computer  program.  The  description  and  use  of  a  direct  shear 
device  for  the  measurement  of  pertinent  soil  properties  are  documented  in 
Appendix  C,  while  Appendix  D  is  a  notation. 
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CHAPTiR  2 


TERRAIN-VEHICLE  INTERACTION  MODEL 


2.1  INTRODUCTION 

The  basic  concepts  of  the  theory  of  terrain-vehicle  interaction  were 
developed  during  the  1950's  by  Bekker  (Reference  4).  By  assuming  various  load 
distributions  along  the  tracks,  Bekker  was  able  to  develop  several  mathemati¬ 
cal  expressions  relating  the  characteristics  of  the  vehicle  and  the  tractive 
effort  of  the  terrain  during  steering.  By  considering  the  lateral  and 
longitudinal  coefficients  of  friction  between  the  track  and  the  ground, 

Hayashi  (Reference  5)  developed  simple  equations  for  practical  analysis  of 
steering  of  tracked  vehicles.  Hayashi's  work,  however,  did  not  include  the 
effect  of  the  centrifugal  forces  on  steering  performance  of  the  vehicle. 

Kitano  and  Jyorzaki  (Reference  6)  developed  a  more  comprehensive  model  for 
uniform  turning  motion  including  the  effects  of  centrifugal  forces.  This 
model,  however,  is  based  on  the  assumption  that  ground  pressure  is  concen¬ 
trated  under  each  road  wheel  and  the  terrain-track  interaction  is  simulated  by 
Coulomb-type  friction.  The  model  given  by  Kitano  and  Jyorzaki  was  extended  by 
Kitano  and  Kuma  (Reference  7)  to  include  nonuniform  (transient)  motion,  but 
the  basic  elements  of  the  terrain-track  interaction  part  of  the  model  were 
retained.  Baladi  and  Rohani  (Reference  1)  developed  a  model  for  uniform 
turning  motion  parallel  to  the  development  by  Kitano  and  Jyorzaki  insofar  as 
the  kinematics  of  the  vehicle  are  concerned.  In  contrast  to  the  development 
by  Kitano  and  Jyorzaki  (Reference  6),  however,  this  model  is  based  on  a  more 
comprehensive  3oil  model.  Baladi  and  Rohani  (References  2  and  8)  extended  the 
WES  terrain- vehicle  model  completed  in  1979  to  include  nonuniform  (transient) 
motion  on  level  terrain.  In  addition,  the  WES  soil  model  was  modified  to 
include  a  nonlinear  failure  envelope  describing  the  shearing  strength  of  the 
terrain  material  (Reference  2)  and  a  nonuniform  (transient)  turning  motion  on 
sloping  terrain  (Reference  3). 

2.2  SOIL  MODEL 

2.2.1  Strength  Components 

One  of  the  most  important  properties  of  soil  affecting  traf f icabi 1 ity  is 
in  situ  shear  strength.  It  has  been  found  experimentally  that  the  shear 
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strength  of  purely  cohesive  soils  is  relatively  independent  of  the  confining 
stress,  but  strongly  affected  by  the  time  rate  of  shearing.  On  the  other 
hand,  the  shear  strength  of  purely  frictional  soils  is  found  to  be  relatively 
independent  of  the  time  rate  of  loading,  but  strongly  dependent  on  confining 
pressure.  The  shearing  resistance  of  most  soils,  however,  is  due  to  both  the 
frictional  and  cohesive  components.  The  cohesive  and  frictional  components  of 
strength  are  usually  added  together  in  order  to  obtain  the  total  shear 
strength  of  the  material.  For  static  loading  (very  slow  rate  of  deformation), 
the  shear  failure  envelope  is  defined  by  the  following  nonlinear  equation: 


where 


tm  =  A  -  M  exp  (-No) 


t..  =  the  maximum  shearing  strength  of  the  material 

M 


(2.1) 


o  =  the  normal  stress 

A  =  the  strength  of  the  material  when  o  is  large 
A  -  M  =  C  =  the  strength  of  the  material  or  cohesion  when  o  =  0 
N  =  a  material  constant 

Equation  2.1  is  shown  graphically  in  Figure  2.1. 

For  a  linear  shear  failure  envelope,  the  following  equation  can  be  used 
in  place  of  Equation  2.1: 


tu  =  C  +  o  tan 
n 


(2.2) 


where  <)>  is  the  angle  of  internal  friction  of  the  material.  Equation  2.2  is 
shown  graphically  in  Figure  2.2. 


Equations  2.1  and  2.2  can  be  combined  in  the  following  relation: 


t  =  [A  -  M  exp  (-No)](1  -  I)+(C  «•  o  tan  $)  I 
m 


(2.3) 


where 


1  for  linear  failure  envelope 
I  = 

0  f or  nonlinear  failure  envelope 

As  noted  previously,  the  shear  strength  of  cohesive  soils  increases  with 
the  increasing  rate  of  loading.  For  the  range  of  loading  rates  asssociated 


I 

i 

I 


to 
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with  the  motion  of  tracked  vehicles,  the  contribution  to  cohesive  strength  due 
to  dynamic  loading  can  be  expressed  as  C^Ll  ~  exP  ( ~  A  A)  3  ,  where  Cd  and  A 
are  material  constants  and  A  is  the  time  rate  of  shearing  deformation.  In 
view  of  Equation  2.  3.  the  dynamic  failure  criterion  takes  the  following  form: 


x  =  [A  -  M  exp  (-No) ]( 1 
m 


I)  +  (C  +  o  tan  <fr)  I  +  Cd  [1  -  exp  (— A A ) ]  (2.4) 


When  A  equals  zero,  the  dynamic  failure  criterion  (Equation  2.4)  reduces  to 
the  static  failure  criterion  (Equation  2.3).  Both  are  shown  graphically  in 
Figure  2.1  for  nonlinear  failure  envelopes  and  Figure  2.2  for  linear  failure 
envelopes . 


2.2.2  Shear  Stress-Shear  Deformation  Relation 


Prior  to  failure,  the  shear  stress-shear  deformation  characteristics  of  a 
variety  of  soils  can  be  expressed  by  the  following  mathematical  expression 
(Reference  9): 


t  = 


0  tm  “ 


'«  * 


(2.5) 


The  behavior  of  Equation  2.5  is  shown  graphically  in  Figure  2.3,  in  which  t 
denotes  shearing  stress,  A  is  shearing  deformation,  and  G  is  the  initial 


shear  stiffness  coefficient.  Substituting  t,.  from  Equation  2.  h  into 

M 


Equation  2.5,  the  shear  stress-shear  deformation  relation  for  soil  becomes 


G{[A  -  M  exp  (-No)](1  -  I)*(C  ♦  o  tan  <j>)  I  +  C  [1  -  exp  (-AA)  ])  A 

t  - - - 7 —  (2.6) 

G |  A |  +  [A  -  M  exp  ( -No )  3(1  -  I)  +  (C  ♦  a  tan  <p)  I  +  Cd  [1  -  exp  (-AA)  ] 

For  purely  cohesive  soils,  N  equals  zero  (<p  =  0  if  the  failure  envelope  is 
linear)  and  t  is  only  a  function  of  A  and  A  .  For  cohesionless  or 
granular  soils,  M  equals  A  (or  C  =  0  if  the  failure  envelope  is  linear), 
is  zero,  and  t  is  a  function  of  A  and  a  .  For  mixed  soils  exhibiting 
shearing  resistance  due  to  both  frictional  and  cohesive  components,  t  is 
dependent  on  A  ,  A  ,  and  o  .  The  qualitative  behavior  of  Equation  2.6  for 
these  three  conditions  is  shown  in  Figure  2.4.  It  should  be  pointed  out  that 
Equation  2.6  reduces  to  the  rigid  plastic  soil  model  often  used  in  mobility 


studies  when  an  extremely  large  value  is  specified  for  G  and  A  is  set  to 
zero. 

An  appropriate  test  for  determining  the  numerical  values  of  the  material 
constants  in  Equation  2.6  is  an  in  situ  direct  shear  test.  A  field  direct 
shear  device  has  been  developed  at  WES  for  this  purpose.  A  description  of 
this  device  and  the  method  of  analysis  of  the  data  obtained  from  the  direct 
shear  test  are  documented  in  Appendix  C  (Reference  8). 

2.3  DERIVATION  OF  TERRAIN-VEHICLE  MODEL 

2.3.1  Boundary  Conditions 

The  geometry  of  the  vehicle  and  the  boundary  conditions  of  the  proposed 
model  are  shown  schematically  in  Figure  2.5.  The  XYZ  coordinates  are  the 
local  coordinate  system  of  which  X  is  always  the  longitudinal  axis  of  the 
vehicle  and  Y  is  a  transverse  axis  parallel  to  the  ground.  These  axes 
intersect  at  the  center  of  geometry  of  the  vehicle  0  .  The  Z  axis  is  a 
vertical  axis  passing  through  the  origin  0  .  The  center  of  gravity  of  the 
vehicle  (CG)  lies  on  the  X  axis  and  is  displaced  by  a  distance  Cx  from  the 

origin.  The  numerical  value  of  Cx  is  assumed  to  be  positive  if  CG  is 

displaced  forward  from  the  center  of  geometry  of  the  vehicle.  The  XY 
coordinates  of  the  instantaneous  center  of  rotation  (ICR)  are  P  +  Cx  and  R  , 
respectively,  where  P  is  the  offset.  The  center  of  rotation  and  the  radius 
of  the  trajectory  of  the  CG  are,  respectively,  CR  and  Rq  .  The  height  of 
the  center  of  gravity  measured  from  ground  surface  is  denoted  by  H  .  The 

lengths  of  the  track-ground  contact,  the  track  width,  and  the  tread  of  the 

tracks  are  L  ,  D  ,  and  B  ,  respectively.  As  shown  in  Figure  2.5,  the 
components  of  the  inertial  forces  F^  in  X  and  Y  directions  are, 
respectively,  F^x  and  Fcx  .  The  weight  of  the  vehicle  is  W  . 

2.3.2  Stress  Distribution  along  the  Tracks 

Two  types  of  stress  (i.e.,  normal  and  shear  stresses)  exist  along  the 
track.  As  indicated  in  Figure  2.5,  the  normal  stresses  under  the  outer  and 
inner  tracks  are  denoted  by  Rj(X)  and  R2(X)  ,  respectively.  The  components 
of  the  shear  stress  in  X  and  Y  directions  are,  respectively,  Tj(X)  and 
Q ,  (X )  for  the  outer  track,  and  T2(X)  and  q2(x)  for  the  inner  track. 

These  stresses  are  dependent  on  the  terrain  type,  vehicle  configuration,  and 
speed  and  turning  radius  of  the  vehicle. 
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The  magnitude  of  normal  stresses  R,(X)  and  R2(X)  can  be  determined  in 
terms  of  the  components  of  the  inertial  force,  the  track  tensions,  and  the 
characteristics  of  the  vehicle  by  considering  the  balance  of  vertical  stresses 
and  their  moments  in  Figure  2.5.  Thus, 


R*  (x) 


R2(x)  = 


W  [l  .  r  _  h  FCY  FCX  dL  Nl(x) 

^  [2  +  6xcx  "  b  ~  6hx  —  +  — w - 

W  [1  .  ,  .  h  FCY  Fcx  dL  N*(x) 

^  2  +  6xCx  +  b  IT  "  6hx  —  +  — - 


(2.7) 


(2.8) 


where 


h  =  H/L 
b  =  B/L 
d  -  D/L 


cx  =  CX/L 


x  =  X/L 
y  =  Y/L 
z  =  Z/L 

Nj(x)  and  N2(x)  «  contributions  due  to  track  tension. 

The  components  of  the  shear  stress  in  the  X  and  Y  directions  along 
both  the  outer  and  inner  tracks  can  be  obtained  by  combining  Equations  2.6, 
2.7,  and  2.8.  Thus  (it  is  noted  that  R,  and  R2  replace  the  normal  stress 
0  in  Equation  2.6), 


Tjtx)  - 


r  /,'nri(x)\l  ri<*> 

|a  -  m  expy - - - yj  (1  -  I),  o  «  — - —  tan  +  I  ♦  c^C  1  - 

.  .  f  /-nr  (x)Yl  7  Vx) 

M|«j|*[a  -  «  expl - j  —  I  (1  -  D*  o  *  tan  ♦  I  ♦  c^t  1 


exp  ( -X6  A ) 3 


exp  (- A6t ) 31 


(2.9) 


L 

m  exp 

t  VJ«  - 

L: 

^  __ 

i(x)  i  •  ) 

-J—  tan  ♦]!  *  odt>  -  e*P  (-iSj)}  / 

»lsi 

/-nr  (x)YT 

a  -  m  exp  ^ - - - )|<1  -  I 

)  * 

Tj (x)tan» 

d 

I  ♦  ed[l  -  exp  (-iijtJ 

sin  T.  (2.10) 


where 


i  =  1,2  (1  for  the  ou;.er  track  and  2  for  the  inner  track) 

r,  (x)  -  dL,V  (x)/W 


<Si  -  Ai/L 
6^  =  A^/L 
M  =  GL*/W 


A  =  AL 


a  =  AL  /W 


ra  =  ML  /W 


n  =  NW/L 


c  -  CL  /W 


c  .  -  Cl  /W 
d  d 


The  variables  Y!  and  Y2  in  Equations  2.9  and  2.10  are  the  slip  angles  and 
can  be  written  as 


Y,  =  tan 


Y2  =  tan 


.  X  -  P  -  C, 


i  X  -  P  -  C, 


=  tan 


i  x  -  p  -  c 

—  4  -v 


(2.11) 


_1  X  -  p  -  c 


where 

Si  =  Cx/L 
C2  =  C2/L 
p  =  P/L 

The  parameter  Ci  is  the  distance  between  the  instantaneous  center  of 
rotation  of  the  outer  track  ICj  and  its  axis  of  symmetry,  and  C2  is  the 
distance  between  the  instantaneous  center  of  rotation  of  the  inner  track  IC2 
and  its  axis  of  symmetry  (Figure  2.6). 

In  order  to  use  Equations  2.9  through  2.11,  the  normal  stress  contribu¬ 
tions  due  to  track  tensions  N,(x)  and  N2(x)  ,  the  track  slip  velocities  and 
displacements  (i.e.,  A,  ,  A,  ,  A2  ,  and  A2)  ,  and  the  inertial  forces  FqX 
and  Fcy  have  to  be  determined.  These  factors  are  discussed  in  the  following 
paragraphs . 


v: 

'v* 

vy 


lii 


K 


& 


2.3.3  Normal  Stress  Contribution  Due  to  Track  Tension 


The  effect  of  track  tension  on  the  normal  stress  distribution  is  influ¬ 
enced  considerably  by  the  motion  of  the  vehicle.  At  relatively  low  speed, 
tractive  effort  is  applied  to  the  outer  track,  while  braking  force  is  applied 
to  the  inner  track  (Figure  2.7a).  At  high  speed,  on  the  other  hand,  tractive 
efforts  are  applied  to  both  tracks  (Figure  2.7b). 

The  angles  0  and  0  in  Figure  2.7  are  the  approach  and  departure 

3  U 

angles  of  the  track  envelope,  respectively.  The  forces  T,  and  T2  are  the 
track  tensions  in  the  outer  and  inner  tracks,  respectively.  These  forces  can 
be  obtained  by  integrating  Equation  2.9.  Thus, 


1 

2 

T,  =  dL2 /  T,(x)  dx 

1 

2 


1 

2 

T2  =  dL  /  T2(x )  dx 
"  2 


(2.12) 


The  normal  stress  distributions  are  influenced,  however,  by  the  vertical 
components  of  the  forces  T,  and  T2  ;  namely,  n,  ,  n2  ,  and  n£  .  The 
values  of  n,  ,  n2  ,  and  n'2  are 


n,  =  T,  sin  0d 


(2.13) 


n2 


T2  sin  0  if  2  0 

3 

0  if  S2  <  0 


(2.14) 


nj 


f-T2  sin  0^  if  52  <  0 


if  52  2  0, 
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(2.15) 


With  the  determination  of  the  forces  nj  ,  na  ,  and  n^  ,  the  normal  stress 
contributions  due  to  track  tension  may  be  determined. 

Since  the  tracks  are  assumed  to  be  rigid,  the  normal  stresses  due  to 
track  tension  may  be  distributed  according  to  the  following  equations 
(Figure  2.8): 


and 


Nj (x)  =  ax  +  m 


i.  1  1 

for  J-Iiisi 


2n, 


Nt(x)  =  ax  +  m  +  — j—  (x  +  ^  ~  f)  for  -  4*xS'r-i 
O  . ,  2  L  2 

d£ 


L  2 


2n, 


1  .  i, 


1  l 


1 


N, (x)  =  ax  +  mT - r  (x  +  -  +  r)  for  ^  -  S  x  ^  - 


dS. 


2  L 


2  L 


N, (x)  =  ax  +  m. 


for  J-isxsi-f 


N,(x) 


2nJ 


ax  +  m  +  — j-  (x  *  •= 
d* 


b  f°>--5SxSI-5 


(2.H 


(2.17 


in  which  i.  is  the  distance  between  two  adjacent  wheels,  and  a  ,  mQ  ,  and 
mj  can  be  determined  by  considering  the  equation  of  equilibrium  of  normal 
stresses  and  the  moments  of  these  stresses.  Thus, 


1_ 

2 


I 


(ax  +  m  )  dx  + 
o 


2Pi 

- 2  (X 


dx 


=  0 


(2.1 


J  (ax 


1 

2 


V 


dx 


2 


£ 

L 


(x  - 


d£ 


dx 


i  x 

L  2 

I 

1 

2 


2nJ 


d£ 


(x  +  ^  +  £)  dx 


and 


1 

2 

J 

x 

2 


1  G, 

(2ax  +  mT  +  mMy  +  r  -  x)  dx 
I  o  2  L 


x 

2 


x 

2 


2ns 

£  d£ 
L 


1 

2 

J 

2 


2(n,  +  n'2) 


d£ 


2  1 
2 -  (X  +  ■= 


L^2  +  L~X^dX_0 


Equations  2.18  through  2.20  contain  three  unknowns:  a  ,  mQ  ,  and 
Completing  the  integrations  results  in 


a  =  -!y  (3  -  ^)<n2  -  n* 
dL  L 


n, ) 


m 


o 


m 


m 


1 


(n2  ♦  nj) 


Lf*.*1'*  f  fc  J  ^  -rJq,  fV’IAt  /  *»,  «  »»•■■»  »•,»  —'V_  «■ 


Substitution  of  Equations  2.21  through  2.23  into  Equations  2.16  and  2.17  leads 


Nj(x)  -  — *7  [(3  -  28)  (n2  -  nj  -  n^  x  +  nj  for  B  -  4  S  x  £  X, 
dL 


Nj  (x ) 


(3-28)  (n2  -  nj  -  nt )  + 


-  x  *  (H1)  n'j 


for  ■  xM  -  ^ 


N2(x)  =  1(3  -  28)  (n2 


2n*  /l  -  8 

-  nj  -  nt)  -  -r-  x  ♦  [— n2  +  nj 


for  ^  -  B  S  x  £  ^ 


N2(x)  =  — [(3  -  28)  (n2  -  ni  -  n, )  x  ♦  n2  ♦  nj] 


for  8  - 


N2(x)  =  — —f  (3  -  28)  (n2  -  nj  -  n,)  +  -j 

dL  L  8 


-  x  +  (H^)  n*  +  ns 


for  -  ^  S  x  S  6  -  ^ 


where 


(2.24) 


(2.25) 


6  =  1 

Note  that  Equations  2.14  and  2.15  dictate  that  either  n2  or  n’2  in  Equations 
2.24  and  2.25  is  zero. 


iy 

A* 

CJ* 

w 

•V 

At 


i 


icfc 


m 


i 


0.0 


■v  V 


2.3.1*  Kinematics  of  the  Vehicle 


A  tracked  vehicle  in  transient  motion  is  shown  schematically  in  Figure  2.9. 
The  XYZ  coordinates  are  the  local  coordinate  systems  that  are  fixed  with 
respect  to  the  moving  vehicle  (also  see  Figure  2.5).  The  origin  0  of  this 
coordinate  system  stays,  for  all  time,  at  a  distance  from  the  center  of 

gravity  of  the  vehicle.  The  4*$  coordinate  system  is  fixed  on  level  ground, 
and  its  origin  coincides  with  the  center  of  gravity  at  time  zero.  The  vehicle 
can  maneuver  on  the  4'4>  plane  and  the  displacements  of  the  center  of  gravity 
of  the  vehicle  from  this  reference  frame  are  fCt)  and  $(t)  . 

The  velocities  vx  and  Vy  (relative  to  the  origin  of  the  ¥«> 
coordinate  system)  as  well  as  the  velocities  v^  and  v  are  related  to  the 
instantaneous  velocity  v  of  the  CG  by 


The  side-slip  angle  a  ,  which  is  the  angle  between  the  velocity  vector  v 
and  the  longitudinal  X  axis  of  the  vehicle,  is  related  to  the  velocities 
vx  and  Vy  as 


_  i 

a  =  tan 


da 

vx  ’  dt 


(2.27) 


The  yaw  angle  w  and  the  directional  angle  0  are  related  to  a  as 


dQ  da)  da 

’  dt  _  dt  dt 


(2.28) 


Substitution  of  Equation  2.27  into  Equation  2.28  leads  to 

(dv 

v  _I  -  v 

X  dt  i 


(2.29) 


The  radius  of  curvature  of  the  trajectory  of  the  center  of  gravity  (i.e., 
the  distance  between  CR  and  CG  (Figures  2.6  and  2.10)  is 


The  coordinates  of  the  trajectory  of  the  center  of  gravity  of  the  vehicle  can 


be  written  as 


( t )  -  -  /  v  cos  0  dt 
0 


*(t)  “  f  v  sin  Q  dt 
0 


(2.31 ) 


The  coordinates  of  the  instantaneous  center  of  rotation  (ICR)  of  the  hull 
in  the  XY  systems  (Xj  ,  Yj)  and  the  instantaneous  radius  of  curvature  (Rj) 
are  (Figures  2.6  and  2.10) 


xi  ■  p  *  cx  -  ',/if  ♦  cx 


h  •  R  ■  »x 


(2.32) 


R:  =  VR  +  P 


The  instantaneous  velocities  of  an  arbitrary  point  e  of  the  hull  are  shown 
in  Figure  2.10  and  can  be  written  as 


v  =  v  +  Y  — 
eX  X  dt 


VeY  VY  (X  "  CXJ  dt 


(2.33) 


'  Wx  +  Y 


30  *  t 


«  -  V  3?] 


2.3.5  Track  Slip  Velocity  and  Displacement 


Assume 


that  vg1  (vgl  =  Ai)  is  the  slip  velocity  of  an.  arbitrary  point 


e,  of  the  outer  track  and  vg2  (vg2  =  A2)  is  the  slip  velocity  at  point 
e2  (e,  and  e2  have  the  same  abscissa)  of  the  inner  track  (Figure  2.6). 
The  X  and  Y  components  of  these  velocities  are 


v  =  c  — 
sXI  1  dt 


C  T  ^ 
'il  L  Hf 


For  the  outer  track  (2.34) 


Vi  ■  <X  -  P  -  CX>  ft  ■  L<X  -  V  3?  - 


I 

& 


•fV 


& 


et 


i 
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For  the  inner  track 


(2.35) 


VsY2  =  VsY 1  ) 

As  indicated  in  Figure  2.11,  the  angular  velocity  du/dt  and  the  value  of  R 
can  be  written  as 


du) 

dt 

R 


—  fv  “V  “V  +  v  ) 
bL  1  XI  sXI  X2  sX2; 


-  (v  “V  +V  “V  ) 

du>  XI  sXI  X2  sX2; 

at 


where 

vx  »  the  velocity  of  the  outer  track  in  X  direction 

v^^  =  the  velocity  of  the  inner  track  in  X  direction 

The  ratio  of  vxi  and  vX2  is  defined  as  the  steering  ratio  e 


(2.36) 


Thus, 


t 


VX1/VX2 


(2.37) 


Substitution  of  Equations  2.32  and  2.37  into  Equation  2.36  leads  to 


V  =  V 

3X1  X2 

l  bL 

-  lvx +  r 

dw\ 

dt/ 

For 

the 

outer 

track 

V  =  V 

sX2  X2 

/  bL 

(  X  2 

duA 

dt/ 

For 

the 

inner 

track 

(2.38) 


Comparison  between  Equations  2.37  and  2.38  and  Equations  2.3*4  and  2.35  results 
in 


5i  - 

£2  = 


1 


(2.39) 


The  slip  velocities  and  displacements  of  the  outer  and  inner  tracks  can  be 
obtained  from  Equations  2.3*4,  2.35,  and  2.38.  Thus, 
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where 


tj  =  (L/2  -  X)/vx1 
t2  =  (L/2  -  X)/vX2 

AI1  =  initial  displacement  of  the  outer  track 

Ai2  =  initial  displacement  of  the  inner  track 

The  values  of  A^  and  Aj2  depend  on  the  balance  between  all  forces  and 
moments  applied  on  the  vehicle  at  zero  velocity.  The  forces  applied  on  the 
vehicle  at  zero  velocity  are  in  turn  dependent  on  the  rolling  resistance. 
Within  the  framework  of  the  present  model,  the  balance  of  forces  and  moments 
dictates  that  the  initial  displacements  be  numerically  equal  to  the  coeffi- 

^11  ^12 

cient  of  rolling  resistance  (i.e.,  — r—  »  — j—  =  ).  The  coefficient  of 

L  lj 

rolling  resistance  must  be  measured  experimentally  or  calculated  from 
empirical  relations  presented  in  Section  2.3.7  (Reference  10). 

2.3.6  Inertial  Forces 

According  to  Figure  2.9,  the  relationship  between  the  velocities  v^  a 

V.  and  the  velocities  v„  and  v 
<P  AY 


can  be  written  as 


The  acceleration  in  ¥  and  ♦  direction,  a^  and  a^  ,  can  be  written  as 


.  .fh 

f  dt 


dv 
J4>  dt 


4> 


(2.-43) 


The  forward  and  lateral  accelerations,  a^  and  a^  ,  can  be  written  in  terms 


of  a^  and  a^  as 


av  =  -a,„  cos  w  +  a^  sin  w 
X  r  9 


a„  =  -a...  sin  <*>  -  a.  cos  w 

It  w 


(2.44) 


Substitution  of  Equations  2.42  and  2.43  into  Equation  2.44  leads  to 


dVX  du> 

ax  =  dt  *  VY  dt 


da 

aY  dt  VX  dt 


(2.4?) 


Hence,  the  X  and  Y  components  of  the  inertial  force  can  be  written  as 

F  -  £  a 
CX  '  g  ax 

P  w 

FCY  g  Y 

where  g  is  the  acceleration  of  gravity. 

2.3.7  Rolling  Resistance 

The  rolling  resistance  R  is  a  function  of  terrain  type,  vehicle  speed, 

s 

track  condition,  etc.  Therefore,  rolling  resistance  should  be  measured  for 


‘'v  b  n 


every  specific  condition.  In  this  formulation,  however,  the  rolling  resis¬ 
tance  is  assumed  to  be  proportional  to  normal  load.  Thus, 


— ~T  1$  /  [iMx)  +  r2(x)]dx 

dL  1 


(2.47) 


As  stated  above,  the  coefficient  of  rolling  resistance  can  be  measured 
experimentally  or  calculated  from  an  empirical  procedure.  Rula  and  Nuttall 
(Reference  10)  presented  such  an  empirical  procedure  by  which  the  coefficient 
of  rolling  resistance  is  calculated  in  terms  of  the  vehicle  characteristics 
and  the  WES  cone  index.  The  procedure  involves  the  following  steps. 

1.  Determine  the  mobility  index  (MI)  for  the  tracked  vehicle  of  interest 
using  the  following  expression: 


■contact  -|  r 
pressure  weight 
■factor  J  [.factor, 
'track  I  Igrouserl 
_factorj  [factor  J 


bogie 

factor 


clear¬ 

ance 

factor 


engine 

factor 


transmi s- 

sion 

factor 


(2.48) 


where 


Contact 
pressure  = 


gross  weight,  lb 


factor  area  of  tracks  in  contact  with  ground,  in. 


(2.49) 


Weight  factor: 


Track  factor  = 


Less  than  50,000  lb  =  1.0 
50,000  to  69,999  lb  =  1.2 
70,000  to  99,999  lb  =  1.4 
100,000  lb  or  greater  =  1.8 


track  width,  in. 


Grouser  factor:  Grousers  less  than  1.5  in.  high  =  1.0 
Grousers  more  than  1.5  in.  high  =  1.1 


A  table  of  factors  for  converting  non-SI  units  of  measurement  to  SI 
(metric)  units  is  presented  on  page  4. 


(2.50) 


(2.51  ) 


(2.52) 


DEKL 


*  ij*  ■  5  * 


(2.53 


Bogie  factor 


_ gross  weight,  lb,  divided  by  10 _ 

(total  number  of  bogies  on  tracks  in  contact 
with  ground)  (area,  in.2,  of  1  track  shoe) 


Clearance  factor 


clearance,  in. 
10 


(2.54 


Engine  factor: 


2  10  hp/ton  of  vehicle  wt  =  1.00 
<  10  hp/ton  of  vehicle  wt  =  1.05 


(2.55 


Transmission 

factor:  Automatic  =  1.0;  manual  =  1.05 


(2.56, 


2.  Determine  the  vehicle  cone  index  VCI1  for  one-pass  traffic  using 
the  expression 


VCI,  -  7.0  .  0.2  HI  - 


(2.57) 


3.  The  coefficient  of  rolling  resistance  is  then  determined  by  the 
following  equation:* 


«-0-0*5*  a ---2;c?!5. 6.5  <2-58> 

where  Cl  is  the  WES  cone  index  for  the  particular  terrain  of  interest.  Note 
that  Cl  must  be  equal  to  or  greater  than  VCI,  in  order  for  the  vehicle  to 
complete  one  pass. 

The  value  of  Cl  must  be  determined  experimentally .  However,  if  such 
measurement  is  not  available,  Cl  can  be  estimated  from  the  parameters  C 
and  <(>  in  the  soil  model.  The  following  empirical  relation  is  often  used  to 
relate  Cl  to  C  and  $  : 

Cl  =  1 2C  (in  psi)  +  4$  (in  degrees)  (2.59) 


In  Reference  10,  the  rating  cone  index  RCI  rather  than  cone  index  Cl 
is  used  to  calculate  .  RCI  is  the  product  of  measured  cone  index  and 
remolding  index  RI  ,  and  is  a  valid  description  only  for  fine-grained  soils 
and  for  sands  with  fines,  poorly  drained.  RI  is  a  ratio  that  expresses  the 
change  in  strength  of  a  fine-grained  soil  or  a  sand  with  fines,  poorly 
drained,  that  may  occur  under  traffic  of  a  vehicle. 


Rohani  and  Baladi  (Reference  11)  developed  more  accurate  relationships  between 
Cl  and  the  mechanical  properties  of  soil  (such  as  C,  $  ,  the  shear 
modulus  G  ,  and  the  density  Y  ).  These  relationships  are 


Cl  =  -C  cot 


+  2  tan  g  (1  +  sin  <j>)Gm 

(|r)  ta"  ♦ 


tan  g  +  tan 
3  -  sin  $ 


8 


where 


[C  ♦  Y(Z 


L)  tan  a  -  [C  ♦  y(2  *  L)  tan  a  »  (2  -  m)YL  tan  »](C  ♦  YZ  tan  »)2 


(2  -  n)  (3  -  m) 


For  granular  materials  C  =  0  and  Equation  2.60  reduces  to 


Cl  = 


2  tan  a(1  +  sin  d> ) G  [  3(  tan 


(|  y)  tan  4> 


3(tan  q  +  tan  4>) 
3  -  sin  4) 


where 


.  £Xil  *  ^  tan  *1^  m  -  [Y  (Z  *  L)  tan  »  *  (2  -  m)  YL  tan  »](YZ  tan  $)2  m 

(2  -  m)  (3  -  in) 


For  cohesive  soil  where  tan  <t>  =  0  ,  the  expression  for  Cl  is  given  by 


Cl  «  4  C(1  ♦  In  5)  +  2k  C  ♦  Y ( 2  ♦  L/3) 
J  L  D 


where 

C 

= 

cohesion 

♦ 

= 

angle  of  internal  friction 

G 

= 

shear  modulus  of  the  material 

Y 

= 

density  of  the  material 

Cl 

= 

cone  index 

D 

s 

diameter  of  the  cone 

L  =  length  of  the  cone 
2a  =  apex  angle  of  the  cone 


(2.60) 


(2.61  ) 


(2.62) 


depth  of  penetration 
(4  sin  $)/[sin  <j>  +  1 )  ] 


m 


In  Equation  2.60  through  2.62,  Z  =  0  corresponds  to  a  fully  embedded  cone. 
These  equations  express  the  cone  index  Cl  in  terms  of  the  soil  parameters 
C  ,  $  ,  Y  ,  and  G  ,  the  depth  of  penetration  Z  ,  and  the  geometry  of  the 

cone. 


2.3*8  Equation  of  Motion 

Steerability  and  stability  of  tracked  vehicles  depend  on  the  dynamic 
balance  between  all  forces  and  moments  applied  on  the  vehicle.  According  to 
Figure  2.5,  the  following  three  equations  govern  the  motion  of  the  vehicle: 


^  l 

2  2 

/  [tt(x)  +  t2(x)]  dx  -  ^  J  [r , (x)  +  r2(x)]dx  =  fcx  (2.63) 

2  2 


2 

/  [qt (x)  +  q2(x)]  dx  *  fCY 
2 


(2.64) 
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-  /  [q,(x)  +  q2(x)]  (x  -  c  )  dx  +  -  /  [t,  (x)  -  t2(x)]  dx 


1_ 

2 


1 
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(2.65) 
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I  w* 
z  d  u> 


2  6  /  Cr2(x)  ~  r,(x)]  dx  =  — — 2 

l_  dt 

2 


where 


.  ,  .  dL 

t,(x)  =  _  Tl(x) 


t2(x)  =  ^  T2(x) 
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r> 


Qi(x)  *  nr-  Q» (x) 


q2(x)  -  %-  Q2(x) 


CX  “  W 


CY  W 

I  =  mass  moment  of  inertia  about  an  axis  passing  through  the  CG  of 
z 

the  vehicle  and  parallel  to  the  Z  axis  (Figure  2.5 ) 

Equations  2.63  through  2.65  are  the  equations  of  motion  for  the  nonuni¬ 
form  (transient)  motion.  As  will  be  shown  below,  the  uniform  turning  motion 
can  also  be  obtained  from  these  equations. 

2.3.9  Treatment  of  Sloping  Terrain 

Figure  2.12  shows  schematically  a  tracked  vehicle  under  nonuniform  (tran¬ 
sient)  turning  motion  on  a  terrain  with  slope  angle  n  .  In  this  case,  the 
weight  of  the  vehicle  W  could  be  resolved  into  a  normal  component  (normal 
to  the  terrain)  WN  and  a  parallel  component  .  Thus, 


=  W  cos  n 


WT  =  W  sin  n 


(2.66) 


In  general,  the  longitudinal  axis  of  the  vehicle  X  makes  an  angle  x  with 
the  component  W^,  (Figure  2.12).  Therefore,  the  component  WT  could  be 
resolved  into  two  components.  The  first  component  WTX  is  parallel  to  the 
X  axis  of  the  vehicle  and  the  second  component  WTy  is  parallel  to  the  Y 
axis.  Thus, 


WTX  =  WT  cos  x  =  w  sin  n  cos  x 
WTY  =  WT  Sln  x  =  W  sin  n  sin  X 


(2.67) 


to 


Eg 


The  angle  x  is  related  to  the  yaw  angle  through  the  following  relation 


X  -  o)  ♦  v 


(2.68) 


i 


where  v  is  a  constant.  The  numerical  value  of  v  depends  on  the  initial 
position  of  the  vehicle  (  v  »  0,  90,  180,  and  270  degrees  corresponds, 
respectively,  to  the  initial  position  of  the  vehicle  at  points  1,  2,  3,  and  4 
in  Figure  2.12).  Substitution  of  Equation  2.68  into  Equation  2.67  leads  to 


W,^  =  W  sin  n  cos  (a  +  v) 


(2.69) 


=  W  sin  n  sin  (u>  +  v) 


In  view  of  Equations  2.66  and  2.69,  the  normal  stresses  under  the  outer  and 
inner  tracks  (Equations  2.7  and  2.8)  become 


_  ,  .  W  cos  n  , 

R,(x)  =  — 2  — 5 — ■  +  6xc  cos 

dL  x 


h  fFCY  .  .  .1 

n  ~  b  hw —  Sin  n  sin^a)  V'J 


(2.70) 


'  QX 

-  6hx  [-jj-  +  sin  n  cos(u  +  v)  +  - 


DLN, (x) 


n  /  \  W  cos  n  z- 

R2(x)  =  j  - — 1  +  6XC  COS 

dL  d  X 


h  f  fcy 

n  +  b  lir 


sin  n  sin(w  +  v) 


(2.71  ) 


-  6hx  +  sin  n  cos(uj  +  v)  +  - 


DLN,(x) 


Equations  2.70  and  2.71  can  be  combined  with  Equations  2.9  through  2.62  to 
develop  the  equations  of  motion  for  a  sloping  terrain.  Thus, 


iv-vv-fv .*oV' •: .a 


2  2 

2  2 

/  [t»(x)  +  ta(x)]  dx  -  /  [r,(x)  +  r2 (x) ]  dx 

2  l 

2  2 

(2.72) 


sin  n  cos  (w  +  v) 
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/  [qi(x)  +  q2(x)]  dx  =  fCY  -  sin  n  sin  (u>  +  v)  (2.73) 

2 
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/  Cq, (x)  +  q2 (x ) ]  (x 

2 

2 


1 
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-  c  )  dx  ♦  /  [t2(x)  -  tj(x)]  dx 


2 
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h  2 

+  ?  h  f  0,(x)  -  r2 (x )]  dx 

2 

2 


\  d^W 
LW  dt2 


(2.74) 


2.3.10  Sprocket  Power 

The  steering  performance  of  a  tracked  vehicle  may  be  limited  either  by 
its  stability  or  by  the  power  available  at  the  sprockets.  The  powers  that 
must  be  available  at  the  inner  and  outer  track  sprockets,  PT1  and  PT2  , 
respectively,  are 


PT1 

— —  «  v 
w  VUg  xl 


2 

2 

/  T, (x)  dx 

2 

2 


(2.75) 
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(2.76) 


=  V  -  /  T2 (x )  dx 


^TEg  -  vx2  J 


where  Tj(x)  and  T2(x)  are  given  by  Equation  2.9.  The  total  power  PT  and 
the  differential  power  PTD  required  at  the  sprockets  are: 


PT  =  PT1  +  PT2 


PTD  =  PT 1  -  PT 2 


(2.77) 


(2.78) 


2.3-11  Power  Output  from  the  Engine 

The  power  that  must  be  available  from  the  engine  is  related  to  the  power 
available  at  the  inner  and  outer  track  sprockets  by  the  following  relation: 


1  PT1  / 
*  2n  e  (£ 


e  +  1  + 


+  PT2  e  +  1  - 


nd  ( e  -  1 ) 


(2.79) 


where 


e  =  steering  ratio 
nd  =  differential  efficiency 

n  =  overall  efficiency  of  the  planetaries  of  the  final  drives 

n  =  pump/ motion  efficiency 
pm 

PT1  =  power  required  by  the  outer  track  sprocket  (Equation  2.75) 

PT2  =  power  required  by  the  inner  track  sprocket  (Equation  2.76) 

The  parameters  ,  q  ,  and  are  vehicle  dependent.  If  the  values 

of  these  parameters  are  not  available,  however,  it  is  recommended  to  use 

n.  =  0.95  ,  n  =  0.95  ,  and  n  =  0.75  . 
d  pm 

2.3*12  Uniform  Turning  Motion 

The  uniform  turning  motion  can  be  obtained  from  the  above  nonuniform 
(transient)  motion  by  setting 


dw  de  .  .  dw 

—  =  —  =  constant;  — , 


(2.80) 


0 


(2.81) 


dv 


dv 


dt 


dt 


Equations  of  motion  for  either  the  nonuniform  or  the  uniform  motions  consti¬ 
tute  three  equations  that  involve  three  unknowns.  The  three  unknowns  are 
either  v^  ,  vy  »  anc*  dw/dt  or  ,  £2  ,  and  p  .  In  order  to  obtain  a 
complete  solution  for  either  of  the  two  sets  of  unknowns,  one  of  the  following 
driving  conditions  must  be  specified: 

1.  Uniform  (steady-state)  motion 

a.  Steering  ratio  for  various  vehicle  velocities. 

b.  Turning  radius  for  various  vehicle  velocities. 

c.  Vehicle  velocity  for  various  steering  ratios. 

d.  Vehicle  velocity  for  various  turning  radii. 

2.  Nonuniform  (transient)  motion 

a.  Time  history  of  the  steering  ratio  e(t)  and  time  history  of  the 
velocity  v(t)  . 

b.  Trajectory  (given  in  terms  of  a  time  history  of  the  instantaneous 
radius  of  curvature  RI(t)  ). 

c.  Time  history  of  the  steering  ratio  e(t)  and  the  velocity  of  the 
vehicle  at  the  beginning  of  the  trajectory. 

d.  Time  history  of  the  velocity  of  the  individual  tracks 
Vxi(t)  and  VX2(t)  . 


A  computer  program  called  TVSTEER  was  developed  to  solve  Equations  2.72 
through  2.7-4  for  both  uniform  (steady-state)  and  nonuniform  (transient) 
motion.  In  addition,  the  program  TVSTEER  contains  an  option  for  calculating 
the  rolling  resistance  (Equations  2.48  through  2.58)  and  the  cone  index 
(Equations  2.60  through  2.62).  A  user's  guide  for  the  program  TVSTEER  is  the 
subject  of  the  remainder  of  this  report. 


SHEARING  OEFORMATION,  A 
a.  PURELY  COHESIVE  SOIL  (r  INDEPENDENT  OF  cr) 


SHEARING  OEFORMATION,  A 
b.  GRANULAR  SOIL  <r  INDEPENDENT  OF  A) 


SHEARING  DEFORMATION,  A 
c  MIXED  SOIL  (r  DEPENDENT  ON  BOTH  a  AND  A) 


Figure  2. A.  Qualitative  behavior  of  the  soil  model 
for  various  types  of  soil. 
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b.  HIGH  SPEED 

Effect  of  track  tension  on  normal 
stress  distribution. 
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2.12.  Effect  of  sloping  terrain  on  transient  motion. 
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USER’S  GUIDE  FOR  TVSTEER 

3.1  INTRODUCTION 

The  computer  program  TVSTEER  is  used  for  predicting  the  steering 
performance  of  high-mobility/agility  tracked  vehicles  in  uniform  turning 
motion  (steady-state)  and/or  nonuniform  (transient)  motion.  It  incorporates 
the  theory  presented  in  Chapter  2  of  this  report  and  uses  a  Newtonian 
iterative  technique  to  solve  the  equations  of  motion.  The  purpose  of  this 
user's  guide  is  to  explain  input  requirements  for  solving  a  particular 
problem,  how  to  run  the  various  options  of  the  program  (i.e.,  steady-3tate, 
transient,  rolling  resistance  calculation,  cone  index  calculation,  etc.),  and 
to  describe  the  output  generated  by  the  program. 

TVSTEER  consists  of  (a)  the  main  program  (see  flowchart  in  Figure  3.1), 
(b)  two  driver  subroutines,  one  called  STEADY  (see  flowchart  in  Figure  3.2) 
for  uniform  turning  motion  and  the  other  called  TRANS  (see  flowchart  in 
Figure  3-3)  for  nonuniform  motion,  (c)  subroutine  AGIL  which  solves  the 
equations  of  motion,  (d)  subroutine  MIVCI  for  computing  the  one-pass  vehicle 
cone  index  (VCI),  and  (e)  subroutine  CONE  for  computing  the  WES  cone  index 
(Cl).  A  complete  program  listing  is  presented  in  Appendix  A.  The  subroutines 
MIVCI  and  CONE  compute  the  VCI  and  Cl  whenever  these  two  quantities  are  not 
specified  in  the  input. 

The  driver  subroutine  STEADY  can  be  utilized  for  each  of  the  following 
input  conditions: 

1.  Input  steering  ratio  and  velocity  increment. 

2.  Input  turning  radius  and  velocity  increment. 

3.  Input  vehicle  velocity  for  various  steering  ratios. 

4.  Input  vehicle  velocity  for  various  turning  radii. 

For  cases  1  and  2,  the  subroutine  STEADY  increases  the  velocity  until  the 
vehicle  becomes  unstable.  In  addition,  if  the  available  power  is  specified, 
the  program  will  produce  and  print  a  complete  set  of  output  when  this  power  is 
reached.  Note  that  in  case  the  subroutine  STEADY  does  not  converge  at  some 
velocity,  it  will  search  for  a  smaller  velocity  increment  at  which  it  can 
converge  (the  smallest  velocity  increment  specified  in  the  code  is  0.03  mph). 
For  case  3,  the  steering  ratio  is  incremented  while  the  vehicle  velocity  is 
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held  constant.  For  case  4 ,  however,  the  turning  radius  is  decremented  while 
the  vehicle  velocity  is  held  constant. 

Similarly,  the  driver  subroutine  TRANS  can  be  utilized  for  each  of  the 
following  driving  conditions: 

1.  Input  the  time  histories  of  the  steering  ratio  and  the  velocity. 

2.  Input  the  time  histories  of  the  turning  radius  and  the  velocity. 

(The  program  can  be  modified,  however,  to  input  the  trajectory  instead  of  the 
time  history  of  the  turning  radius.) 

3.  Input  the  time  history  of  the  steering  ratio  and  the  velocity  of  the 
vehicle  at  the  beginning  of  the  trajectory. 

H.  Input  the  time  histories  of  the  individual  track  velocities. 

In  cases  1  and  2,  the  actual  velocity  used  by  the  program  could  be  less  than 
the  input  velocity  if  the  calculated  power  exceeds  a  specified  available 
power.  In  case  3.  the  velocity  of  the  vehicle  is  adjusted  to  keep  the 
required  power  between  95  and  100  percent  of  the  available  power.  If  the 
vehicle  becomes  unstable  during  this  process,  the  velocity  and/or  the  steering 
ratio  will  be  adjusted  in  an  attempt  to  stabilize  it.  Note  that  all  time 
history  input  can  be  described  by  mathematical  functions  or  can  be  digitized. 

This  user's  guide  contains  a  description  of  typical  input  and  output,  a 
glossary  of  the  important  variables,  flowcharts  (Figures  3.1 “3. 3),  a  listing 
of  the  program  TVSTEER  (Appendix  A),  and  several  sample  runs.  Figures  3. 4  and 
3.5  illustrate  sample  inputs  for  a  steady-state  run  and  a  transient  run, 
respectively.  A  total  of  ten  sample  output  are  illustrated  in  Figures  3.6- 
3.15.  Figures  3.6  through  3.9  contain  sample  output  for  a  steady-state  motion 
corresponding  to  cases  1  through  4 ,  respectively.  Figures  3.10  through  3. 13 
portray  sample  output  for  a  transient  motion  corresponding  to  cases  1  through 
4 ,  respectively.  All  the  above  examples  are  calculated  using  a  linear  soil 
model  and  zero  track  tension.  To  demonstrate  the  effect  of  the  track  tension, 
case  1  of  the  steady-state  motion  was  repeated  and  the  results  are  illustrated 
in  Figure  3-14.  Also,  case  1  of  the  steady  state  motion  was  repeated  for  a 
nonlinear  soil  model.  The  result  of  this  case  is  shown  in  Figure  3.15. 

Program  TVSTEER  has  been  coded  in  Honeywell  FORTRAN  77  for  the  time¬ 
sharing  subsystem  of  the  Honeywell  DPS-8  digital  computer  currently  in 
operation  at  WES.  The  processor  time  for  each  sample  run  is  shown  at  the  end 
of  the  output. 
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3.2  INPUT 

TVSTEER  requires  input  from  the  keyboard  describing  the  vehicle,  the  soil 
properties  along  with  the  slope  of  the  terrain  and  the  directional  angle  of 
the  vehicle,  and  the  driving  conditions  (i.e.,  steady-state  or  transient 
motions,  cases,  etc.).  Note  that  the  program  could  be  modified  to  read  input 
from  a  file.  All  input  except  the  radius  of  curvature  and  the  velocity  are  in 
the  non-SI  units  of  inches,  seconds,  pounds,  and  degrees.  The  input  for  the 
radius  of  curvature  and  the  velocity  are  in  feet  and  miles  per  hour,  respec¬ 
tively.  These  quantities,  however,  are  converted  internally  in  the  program  to 
inches  and  inches  per  second. 


3.2.1  Input  to  the  Main  Program 

The  first  line  of  input  contains  the  test  title: 

TESTN  —  Identification  label  of  up  to  40  alphanumeric  characters 
The  second  line  of  input  contains  the  following  vehicle  parameters: 

WT  =  Weight  of  the  vehicle,  lb. 

L  =  Length  of  the  track  in  contact  with  the  ground,  in. 

H  -  Height  of  the  center  of  gravity,  in. 


& 


D  M 

B  = 

SL  = 

CX  = 
THETAD  = 

THETAA  = 

IZ  - 

VCI1  - 


Track  width,  in. 

Track  tread  (distance  between  the  center  lines  of  the 
tracks),  in. 

Distance  between  the  center  lines  of  adjacent  wheels  (if 


the  model  is  to  be  run  without  using  track  tension  set 
SL  =  0),  in. 

Abscissa  of  the  center  of  gravity  of  the  vehicle,  in. 
Departure  angle  of  the  track  envelope  (THETAD  is  used 
only  when  running  with  track  tension),  deg. 

Approach  angle  of  the  track  envelope  (THETAA  is  used  only 
if  running  with  track  tension  and/or  computing  VCI1), 


Mass  moment  of  inertia  of  the  vehicle  about  an  axis 
passing  through  its  center  of  gravity  and  parallel  to  the 

p 

Z  axis  (not  needed  for  steady-state),  in-lb/sec  . 

Vehicle  cone  index  for  one-pass  traffic.  (VCI1  may  be 
computed  using  subroutine  MIVCI  by  inputting  VCI1  »  0. 
MIVCI  will  ask  for  additional  input  as  explained  below.) 


The  third  line  of  input  contains  the  following  variables  which  are  needed 
by  subroutine  MIVCI: 

GH  -  Height  of  the  grouser,  in. 

NB  »  Total  number  of  bogies  on  track  in  contact  with  ground. 

TSL  -  Track  shoe  length,  in. 

HP  -  Specified  engine  horsepower,  hp. 

TRNT  =  Type  of  transmission  (0  for  manual,  1  for  automatic). 

Note  that  the  above  line  of  input  is  skipped  if  VCI1  is  specified. 

The  fourth  line  of  input  contains  the  following  soil  model  parameters: 

Cl  =  WES  cone  index  (Cl  may  be  computed  using  subrouting  CONE 
by  inputting  Cl  •  0,  subroutine  CONE  will  ask  for 
additional  input  as  explained  below). 

A  «  Material  constant  in  failure  envelope  (Equals  C  for 
linear  soil  model),  psi. 

SM  =  Material  constant  in  failure  envelope  (Equals  0  for 
linear  soil  model),  psi. 

SN  =  Material  constant  in  failure  envelope,  1/psi. 

SXI  =  Variable  used  to  determine  the  type  of  soil  model  used  (0 
for  nonlinear  failure  envelope,  1  for  linear  failure 
envelope) . 

CD  *  Added  cohesive  strength  due  to  dynamic  loading,  psi. 

CLAMDA  =  Material  constant  related  to  rate  effect,  sec/in. 

PHI  =  Angle  of  internal  friction,  deg. 

G  =  Shear  modulus,  psi /in. 

SF  =  Coefficient  of  rolling  resistance.  (SF  will  be  computed 
by  the  code  when  SF  =  0.) 

ETA  =  Angle  of  sloping  terrain,  deg. 

CHI  =  Directional  angle  of  the  vehicle  on  sloping  terrain  (for 
transient  motion  CHI  is  the  initial  angle),  deg. 

C  =  Static  cohesive  component  of  shear  strength  (C  is 
computed  by  code  as  C  =  A  -  SM) ,  psi. 

The  fifth  line  of  input  contains  the  following  variables  which  are  needed 
by  subroutine  CONE: 

CL  =  Cone  length,  in. 

DI  =  Cone  diameter,  in. 


GAMA  =  Density,  psi/in. 

2  =  Depth  of  penetration,  in. 

Note  that  the  above  line  of  input  is  skipped  if  Cl  is  specified. 

The  last  line  of  input  for  the  main  program  is  the  type  of  run: 

IRUN  =  Run  type  (0  for  steady-state,  1  for  transient). 

3.2.2  Input  to  Subroutine  STEADY  (Steady-State  Driver) 

The  first  line  of  input  is: 

ICASE  =  Case  number:  1 ,  2,  3.  ^  or  ?  (to  have  an  input  menu 
displayed,  enter  ?). 

The  content  of  the  second  line  of  input  is  determined  by  the  case  number 
as  follows: 

1 .  For  ICASE  =  1  : 

E  =  Steering  ratio 

DV  =  Velocity  increment,  mph  (converted  internally  to  in/sec). 
HPAV  =  Horsepower  available  (if  HPAV  =  0,  the  power  check  is 
skipped),  hp. 

2.  For  ICASE  =  2: 

R0  =  Radius  of  the  trajectory  (turning  radius)  of  the  center 
of  gravity  of  the  vehicle,  ft  (converted  internally  to 
in) . 

DV  =  See  ICASE  =  1. 

HPAV  =  See  ICASE  =  1 . 

3.  For  ICASE  =  3: 

E  =  Initial  steering  ratio. 

DE  =  Steering  ratio  increment. 

NP  =  Number  of  points  that  will  be  calculated. 

4.  For  ICASE  =  4: 

R0  =  Initial  radius  of  the  trajectory  (turning  radius)  of  the 
center  of  gravity  of  the  vehicle,  ft  (converted 
internally  to  in). 

DR  =  Turning  radius  decrement,  ft. 

NP  =  Number  of  points  that  will  be  calculated. 

3.2.3  Input  to  Subroutine  TRANS  (Transient  Driver) 

The  first  line  of  input  contains  the  following  information: 

DT  =  Time  increment,  sec. 


jfrf  ^ 


IPRINT 


Number  of  points  to  be  calculated. 
Print  skip  increment. 


The  second  line  of  input  is  the  case  number: 

ICASE  -  Case  number:  1,  2,  3.  or  ?  (to  have  an  input  menu 
displayed,  input  ?). 

The  remaining  lines  of  input  are  determined  by  case  number: 

1 .  For  ICASE  =  1 : 

The  first  line  of  input  contains  the  following  power  and  velocity 


information: 


V3,  V4 


-  Horsepower  available  (if  HPAV  =  0,  the  power  check  is 
skipped),  hp. 

*  Initial  velocity,  mph. 

=  Maximum  change  in  velocity,  mph. 

*=  Coefficients  which  are  related  to  the  rate  of  change  in 
velocity.  For  a  constant  velocity,  input  V2,  V3,  and  V4 
as  zero. 


The  second  line  of  input  is  the  values  of  the  coefficients  in  the 
equation  of  the  steering  ratio: 

El  »  Initial  steering  ratio. 

E2  =  Maximum  change  in  the  steering  ratio. 

E3,  E4  =  Coefficients  which  are  related  to  the  rate  of  change  in 

the  steering  ratio.  For  a  constant  steering  ratio,  input 
E2,  E3,  and  E4  as  zero. 

2.  For  ICASE  =  2: 

The  first  line  of  input  is  the  same  as  the  first  line  of  input  for 
ICASE  =  1. 

The  last  line  of  input  is  the  values  of  the  coefficients  in  the  equation 
of  the  turning  radius: 

RI1  =  Initial  instantaneous  radius  of  curvature,  ft. 

RI2  =  Maximum  change  in  this  radius,  ft. 

RI3,  RI4  =  Coefficients  which  are  related  to  the  rate  of  change  in 
this  radius.  For  a  constant  radius,  input  RI2,  RI3,  and 
RI4  as  zero. 

3.  For  ICASE  «  3: 

The  only  line  of  input  is: 
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HPAV  -  Horsepower  available,  hp. 

VI  =  Initial  velocity,  mph. 

For  this  case,  the  digitized  time  history  is  stored  in  the  array  EP  in  a  data 
statement. 

1*.  For  ICASE  = 

The  first  line  of  input  is: 

VX11  *  Initial  velocity  of  the  outer  track,  mph. 

VX1 2  =  Maximum  change  in  this  velocity,  mph. 

VX13*  VX1  4  =  Coefficients  which  are  related  to  the  rate  of  change  in 

this  velocity.  For  a  constant  velocity,  input  VX12, 

VX1 3,  and  VX1 4  as  zero. 

The  last  line  of  input  is: 

VX21  =  Initial  velocity  of  the  inner  track,  mph. 

VX22  =  Maximum  change  in  this  velocity,  mph. 

VX23,  VX24  =  Coefficients  which  are  related  to  the  rate  of  change  in 
this  velocity.  For  a  constant  velocity,  input  VX22, 

VX23,  and  VX24  as  zero. 

3.3  OUTPUT 

The  output  produced  by  TVSTEER  is  presented  in  tabulated  form  in  this 
user's  guide  (see  Figures  3-6-3.15).  Because  all  the  output  values  are  stored 
in  arrays,  a  plot  routine  accessing  these  arrays  could  be  added  easily.  Also, 
a  postprocessor  could  be  added  to  convert  all  the  output  quantities  from  non- 
SI  to  SI  units. 

The  first  part  of  the  output  contains  the  type  of  test  and  most  of  the 
input  as  well  as  the  vehicle  and  soil  model  parameters. 

The  remainder  of  the  output  is  a  table  containing  the  following  variables 
(listed  as  they  appear  in  the  table): 

PT  =  Point  number. 

T  =  Time  (not  tabulated  for  steady-state),  sec. 

P  =  Offset,  P/L 
E  =  Steering  ratio. 

R0  =  Radius  of  the  trajectory  (turning  radius)  of  the  center 
of  gravity  of  the  vehicle,  ft. 

=  Velocity  of  the  vehicle,  mph. 


V 


VX1  =  Longitudinal  component  of  velocity  of  the  outer  track, 
mph. 

VX2  =  Longitudinal  component  of  velocity  of  the  inner  track, 
mph. 

VS1  =  Longitudinal  component  of  slip  velocity  of  the  outer 
track,  mph. 

VS2  =  Longitudinal  component  of  slip  velocity  of  the  inner 
track,  mph. 

W  =  Yaw  angle  (not  tabulated  for  steady-state),  deg. 

WD  =  Yaw  rate,  deg/sec. 

FCX1  =  Forward  acceleration,  g's. 

FCY^  =  Lateral  acceleration,  g's. 

PTE  =  Power  required  from  the  engine,  hp. 

PTS  =  Total  power  required  by  the  sprockets,  hp. 

DXT  =  Abscissa  of  trajectory  (not  tabulated  for  steady-state), 
ft. 

DYT  =  Ordinate  of  trajectory  (not  tabulated  for  steady-state), 
ft. 

The  following  variables  are  available  but  are  not  included  in  the 
tabulated  output. 

PT1  =  Power  required  by  the  sprocket  of  the  outer  track,  hp. 

PT2  =  Power  required  by  the  sprocket  of  the  inner  track,  hp. 

VX  =  Longitudinal  component  of  velocity  of  the  vehicle,  mph. 

VY  =  Lateral  component  of  velocity  of  the  vehicle,  mph. 


1 


There  two  quantities  are  also  used  in  chapter  3  as  the  longitudinal  and 
lateral  components  of  inertial  forces,  respectively. 


Figure  3.1.  Flowchart  for  TVSTEER,  main  program. 


INPUT  E,  DE,  NP 
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PRINT  TITLES 
AND  ECHO  INPUT 


Figure  3.3.  (Sheet  3  of  3). 
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INPUT  TEST  NO.  <Uf  TO  <0  CHARACTERS)  -  SAflRLE  RUN  FOR  USER'S  GUIDE 
TO  NAME  CORE  COMPUTE  UC 1 1 .  INRUT  VC  11=0. 

INRUT  U T  >  L  >  H , D , D , SL  >  CX , THE  T  AD . THE T AA  > I Z . VC  1 1  =  18000. . 105. ,35.7, 15. t?0. .0. .0. .0. ,30.  .0.  .0. 


FOR  TRNT,  0  =  MANUAL >  I  =  AUTOMATIC 
INPUT  GH.NB.TSL.KP.TRNT  =  1 . 4 ,5.9. 6. 300. . 1 

TO  HAVE  COPE  COMPUTE  CONE  INDEX.  INPUT  CI=0. 

TC  HAVE  CODE  COMPUTE  ROLLING  RESISTANCE.  INPUT  SF=0. 

INPUT  CI.A.SM.SN.SXI.CO.CLAMDA.PHI.G.SF.ETA.CHI  =  0. . 4. .0. .0. .1.. >0.18.25.1  SO. .0..C..0. 

TO  SIMULATE  STANDARD  UES  CONE.  INPUT  CL=I . 48, Dt=0. 7? 

FOR  TYPICAL  GAMA, 2.  INPUT  GAMA=0.07.Z=&. 

INPUT  CL. 01. GAMA. Z  =  1  .-18.0.7?. 0.0’»6. 

input  punttre: 

FOR  STEADY-STATE  RUN.  INPUT  0 
FOR  TRANSIENT  RUN,  INPUT  1 
INPUT  YOUR  CHOICE  =  0 

INPUT  CASE!  1,2, 3, A,  OR  ?  =  ? 

case  selections: 

CASE  C1J . VARY  V  (VELOCITY)  i  INPUT  E  (STEERING  RATIO) 

CASE  C2I . VARY  V  (VELOCITY)  INPUT  RO  (TURNING  RADIUS) 

CASE  C 3) . .  . .  VARY  E  (STEERING  RATIO)  !  FIXED  V  (VELOCITY) 

CASE  CO. ..VARY  RO  (TURNING  RADIUS)  ?  FIXED  V  (VELOCITY) 

INPUT  YOUR  CHOICE.'  =  1 
INPUT  E.DV.HPAV  =  1. Art. .1009. 


Figure  3.4.  Sample  input  for  subroutine  STEADY 
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FRN 

INPUT  TEST  NO.  (LIP  TO  <0  CHARACTERS)  *  SAMPLE  RUN  FOR  USER'S  GUIDE 
TO  HAVE  COPE  CONFUTE  VCI1.  INPUT  VCI1=0. 

INPUT  UT.L.H.P.B.SL.CX.THETAP.THETAA.IZ.VCII  -  18000. .103. .33. ?. IS.  >90. >0. »0. .0. ,30. >92009. >20.9 

TO  HAVE  COPE  COMPUTE  CONE  INDEX.  INPUT  CI*0. 

TO  HAVE  COPE  COMPUTE  ROLLING  RESISTANCE.  INPUT  SF-O. 

INPUT  Cl .A.SH.SN.SXI .CP.CLAMDA.PHI .G.SF.E1 A.CHI  =  I6i . ,4. ,0. .0. , 1 . .4. ,0. 18.23. ,30. ,0. *0.  .0. 

INPUT  RUNTTPE1 

FOR  STEAPT-STATE  RUN.  INPUT  0 
FOR  TRANSIENT  RUN,  INPUT  1 
INPUT  YOUR  CHOICE  *  1 

MAXIMUM  POINTS  NP  SET  TO  100 
INFUT  PT.NP.IPRINT  =  0.2.40.1 

INPUT  CASE!  I. 2.3. 4.  OR  T  *  » 

CASE  selections: 

CASE  C1J . INFUT  V.E  (VELOCITY, STEERING  RATIO) 

CASE  [21 . INPUT  V.RJ  (VELOCI  TY.  TURNING  RADIUS) 

CASE  C3J.. .INPUT  E.VI  (STEERING  RATIO-INITIAL  VELOCITY)  i  VARY  V  (VELOCITY) 
CASE  C43 . INPUT  VXI.VX2  (TRACK  VELOCITIES) 

INPUT  YOUR  CHOICE:  =  1 

INFUT  HPAV.V1.V2.V3.V4  =  1000. .5. . -15. . 1 . . 2. 

INPUT  El.E2.E3.E4  =  1.1, -0.4-1. .2. 


Figure  3.5.  Sample  input  for  subroutine  TRANS. 
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Figure  3.7. 


Sample  run  for  steady-state  motion.  Case  2;  soil  model  with 
linear  shear  failure  envelope;  zero  track  tension. 
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Figure  3.8.  Sample  run  for  steady-state  motion.  Case  3;  soil  model  with 
linear  shear  failure  envelope;  zero  track  tension. 
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Figure  3.9. 


Sample  run  for  steady-state  motion,  Case  4;  soil  model  with 
linear  shear  failure  envelope;  zero  track  tension. 
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Figure  3.10.  Sample  run  for  transient  motion,  Case  1;  soil  model  with 
linear  shear  failure  envelope;  zero  track  tension. 
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Figure  3.11.  Sample  run  for  transient  motion,  Case  2;  soil  model  with 
linear  failure  envelope;  zero  track  tension. 
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Figure  3.12.  Sample  run  for  transient  motion.  Case  3;  soil  mod*.:.  i 
linear  shear  failure  envelope;  zero  track  tension. 
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Figure  3.13.  Sample  run  for  transient  motion.  Case  4;  soil  model  with 
linear  shear  failure  envelope;  zero  track  tension. 
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Figure  3.15.  Sample  run  for  steady-state  motion.  Case  1;  soil  model  with 
nonlinear  shear  failure  envelope;  zero  track  tension. 
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APPENDIX  A:  PROGRAM  LISTING 


*m*m*mm*mmmmmmm*m**mmmm*m* 


TVSTEER  IS  A  MODEL  USED  FOR  PREDICTING  STEERING 
PERFORMANCE  OF  HIGH-MOBILITY/AGILITY  TRACKED  VEHICLES 
IN  STEADY-STATE  TURNING  MOTION  AND/OR  TRANSIENT  MOTION. 


INPUT  FOR  TVSTEER 


TESTN  —  TEST  TITLE  (40  CHARACTERS  MAX) 

IRUN  —  RUNTYPE  (0  =  STEADY-STATE*  1  =  TRANSIENT) 


VEHICLE  PARAMETERS 


UT 

L 

H 

D 

B 

SL 

CX 

THETAD 

THETAA 

IZ 


VEHICLE  WEIGHT  (LB) 

CONTACT  LENGTH  OF  TRACK  (IN) 

HEIGHT  OF  THE  CENTER  OF  GRAVITY  (IN) 

TRACK  WIDTH  (IN) 

TRACK  TREAD*  DISTANCE  BETWEEN  CENTER  LINES  (IN) 
DISTANCE  BETWEEN  CENTERS  OF  ADJACENT  WHEELS  (IN) 
ABSCISSA  FOR  CENTER  OF  GRAVITY  (IN) 

DEPARTURE  ANGLE  OF  TRACK  ENVELOPE  (DEG) 

APPROACH  ANGLE  OF  TRACK  ENVELOPE  (DEG) 

MASS  MOMENT  OF  INERTIA  (IN  LB/SEC**2>> 

NOT  APPLICABLE  FOR  STEADY  STATE 
ONE  PASS  VEHICLE  CONE  INDEX 
TO  HAVE  CODE  COMPUTE  VCI1* 

THE  FOLLOWING  INPUT  IS  REQUIRED  BY  SUBROUTINE  MIVCI 


GROUSER  HEIGHT  (IN) 

TOTAL  NUMBER  OF  BOGIES  ON  TRACK 
TRACK  SHOE  LENGTH  (IN) 

HORSEPOWER 

TRANSMISSION  TYPE  ( 0=MANUAL *  1=AUT0MATIC) 


SOIL  MODEL  PARAMETERS 


WES  CONE  INDEX 

TO  HAVE  CODE  COMPUTE  Cl* 

THE  FOLLOWING  INPUT  IS  REQUIRED  BY  SUBROUTINE  CONE 


CL 

DI 

GAMA  - 
Z 


-  CONE  LENGTH  (IN) 

-  CONE  DIAMETER  (IN) 

-  DENSITY  (PSI/IN) 

-  DEPTH  (IN) 


TRACKOOl 

TRACK002 

TRACK003 

TRACK004 

TRACKOOS 

TRACK006 

TRACK007 

TRACK008 

TRACK009 

TRACK010 

TRACK011 

TRACK012 

TRACK013 

TRACK014 

TRACK015 

TRACK016 

TRACK017 

TRACK018 

TRACK019 

TRACK020 

TRACK021 

TRACK022 

TRACK023 

TRACK024 

TRACK025 

TRACK026 

TRACK027 

TRACK028 

TRACK029 

TRACK030 

TRACK031 

TRACK032 

TRACK033 

TRACK034 

TRACK035 

TRACK036 

TRACK037 

TRACK038 

TRACK039 

TRACK040 

TRACK041 

TRACK042 

TRACK043 

TRACK044 

TRACK045 

TRACK046 

TRACK047 

TRACK048 

TRACK049 

TRACK050 

TRACK051 


UWUl 


5? 


5 

‘5 

Is 

s 


A 

SM 

SN 

SXI 

CD 

CLAMDA  — 
PHI 


SHEAR  STRENGTH  OF  SOIL  (PSI) 

MATERIAL  CONSTANT  IN  FAILURE  ENVELOPE  (PSI) 

MATERIAL  CONSTANT  IN  FAILURE  ENVELOPE  (1/PSI) 

TYPE  OF  SOIL  MODEL  (0  =  NONLINEAR.  1  =  LINEAR) 

ADDED  COHESIVE  STRENGTH  DUE  TO  DYNAMIC  LOADING  (PSI) 
MATERIAL  CONSTANT  RELATED  TO  RATE  EFFECT  (SEC/IN) 
ANGLE  OF  INTERNAL  FRICTION  (DEG) 

SHEAR  MODULUS  (PSI/IN) 

ROLLING  RESISTANCE  (MAY  BE  COMPUTED  BY  CODE) 

ANGLE  OF  SLOPING  TERRAIN  (DEG) 

DIRECTIONAL  ANGLE  ON  SLOPING  TERRAIN  (DEG) 

COHESION  (PSI).  C  =  A-SM 


INPUT  FOR  STEADY 


ICASE  --  TYPE  OF  RUN 

E  --  STEERING  RATIO 

DV  --  VELOCITY  INCREMENT  (MPH) 

HPAV  —  AVAILABLE  HORSEPOWER 

RO  --  TURNING  RADIUS  (FT) 

DE  —  STEERING  RATIO  INCREMENT 

NP  --  NUMBER  OF  POINTS 

DR  —  TURNING  RADIUS  INCREMENT  (FT) 


INPUT  FOR  TRANS 


DT  —  TIME  INCREMENT  (SEC) 

NP  —  NUMBER  OF  POINTS 
IFRINT  —  PRINT  SKIP  INCREMENT 
ICASE  --  TYPE  OF  RUN 
HPAV  —  AVAILABLE  HORSEPOWER 
VI  —  INITAIL  VELOCITY  (MPH) 

EP  —  ARRAY  CONTAINING  STEERING  RATIO  TIME  HISTORY 


DRIVER  COEFFICIENTS 

V1-V4.  E1-E4.  RII-RI4.  VX11-VX14.  VX21-VX24 


OUTPUT 


POINT  NUMBER 

TIME  (SEC).  NOT  APPLICABLE  FOR  STEADY-STATE 
OFFSET 

STEERING  RATIO 

RADIUS  OF  TRAJECTORY  AT  CENTER  OF  GRAVITY  (FT) 
VEHICLE  VELOCITY  (MPH) 

OUTER  TRACK  VELOCITY  (MPH) 

INNER  TRACK  VELOCITY  (MPH) 

OUTER  TRACK  SLIP  VELOCITY  (MPH) 

INNER  TRACK  SLIP  VELOCITY  (MPH) 

YAW  ANGLE  (DEG).  NOT  APPLICABLE  FOR  STEADY-STATE 


TRACK052 
TRACK053 
TRACK054 
TRACK055 
TRACK056 
TRACK037 
TRACK058 
TRACK059 
TRACK060 
TRACK061 
TRACK062 
TRACK063 
TRACK064 
TRACK065 
TRACK066 
TRACK067 
TRACK068 
TRACK069 
TRACK070 
TRACK071 
TRACK072 
TRACK073 
TRACK074 
TRACK075 
TRACK076 
TRACK077 
TRACK078 
TRACK079 
TRACK080 
TRACK081 
TRACK082 
TRACK083 
TRACK084 
TRACK085 
TRACK086 
TRACK087 
TRACK088 
TRACK089 
TRACK090 
TRACK091 
TRACK092 
TRACK093 
TRACK094 
TRACK095 
TRACK096 
TRACK097 
TRACK098 
TRACK099 
TRACK100 
TRACK101 
TRACK102 
TRACK  103 
TRACK104 
TRACK105 
TRACK106 
TRACK107 


'MKXMXMi  Ksoaa&e  w max* 


liUSkViJM'i 


LViJ’iia'i 


YAW  RATE  ( DEG/SEC ) 

LONGITUDINAL  ACCELERATION  (G) 

LATERAL  ACCELERATION  <G> 

ROWER  REQUIRED  FROM  ENGINE  (HP) 

POWER  REQUIRED  AT  SPROCKET  (HP) 

ABSCISSA  OF  TRAJECTORY  (FT). 

NOT  APPLICABLE  FOR  STEADY-STATE 
ORDINATE  OF  TRAJECTORY  (FT). 

NOT  APPLICABLE  FOR  STEADY-STATE 
POWER  REQUIRED  BY  SPROCKET  OF  OUTER  TRACK  (HP) 

POWER  REQUIRED  BY  SPROCKET  OF  INNER  TRACK  (HP) 

PT1  AND  PT2  ARE  NOT  TABULATED 
LONGITUDINAL  COMPONENT  OF  VEHICLE  VELOCITY  (MPH) 
LATERAL  COMPONENT  OF  VEHICLE  VELOCITY  (MPH) 

VX  AND  VY  ARE  NOT  TABULATED 


PARAMETER  (  N2  =  100  ) 

REAL  L. 12 

CHARACTER  TESTN*40,DATE*8 

DATA  GRAV.DRAD. VMP  /  386.4,  0.017453293,  0.056818182  / 
DATA  MPR.RDEG  /  0040075040021,  57.29578  / 


COMMON  /ARRAYS/ 


COMMON 

COMMON 

COMMON 

COMMON 

COMMON 


/INOUT  / 
/IN  / 
/VPARAM/ 
/SPARAM/ 
/OUT  / 


COMMON  /MISC 


TP(N2)»VXP(N2) » VYP(N2) ,FCXP(N2) ,FCYP(N2) » 
VSX1P(N2 ) . VSX2P(N2) ,R0P(N2) ,PTEP(N2  >  »PTSP(N2) , 
DXTP(N2),DYTP(N2),WDP(N2),UF(N2) ,VX1P(N2), 
VX2P(N2) ,EP(N2> ,VP(N2) »PP(N2) »PT1P(N2) »PT2P(N2) 
E»VX1 ,VX2,V»RI0L,L0GRE,LVX 
DTIM1 , DTIM2, WPRV, WDPRV, VXPRV, VYPRV, TESTN,  DATE 
WT, L ,H, D , B,3L , CX , THETAA , THETAD, 12 ,  VCI 1 
CI,A,SM,SN,C,SXI,CD,CLAMDA,PHI,G,SF,ETA,CHI 
FCX.FCY, VSX1 , VSX2,PTE,PTS, WD , U, VX, VY, 

P,21 ,22, IDIR, VXD, VYD.PT1 ,PT2 
GRAV.DRAD, VMP, RDEG 


THESE  ARE  CALLS  TO  SYSTEM  DEPENDENT  SUBROUTINES 


CALL  FPARAM(3,MPR) 
CALL  FPARAM( 1 ,132) 
CALL  FTIME(PTI) 

CALL  DATIM(DATE.TIM) 


END  SYSTEM  DEPENDENT  CALLS 


PRINT  *,  'INPUT  TEST  NO.  (UP  TO  40  CHARACTERS)' 

READ  240,  TESTN 
PRINT  *, 

PRINT  *,  'TO  HAVE  CODE  COMPUTE  VCI1,  INPUT  VCI1=0.' 
PRINT  *,  'INPUT  WT,L,H,D,B,SL,CX, THETAD, THETAA, 12, VCI1' 
READ  *,  UT,L,H,D,B,SL,CX, THETAD, THETAA, 12, VCI1 
PRINT  *, 

IF  (VCI1  ,LE.  0.)  CALL  MIVCI 


TRACK108 
TRACK109 
TRACKUO 
TRACK1 1 1 
TRACK112 
TRACK113 
TRACK114 
TRACK115 
TRACK116 
TRACK117 
TRACK118 
TRACK119 
TRACK120 
TRACK121 
TRACK122 
TRACK123 
TRACK124 
TRACK125 
TRACK126 
TRACK127 
TRACK128 
TRACK129 
TRACK130 
TRACK131 
TRACK132 
TRACK133 
TRACK134 
TRACK135 
TRACK136 
TRACK137 
TRACK138 
TRACK139 
TRACK140 
TRACK141 
TRACK142 
TRACK143 
TRACK144 
TRACK145 
TRACK146 
TRACK147 
TRACK148 
TRACK149 
TRACK150 
TRACK151 
TRACK152 
TRACK153 
TRACK154 
TRACK155 
TRACK156 
TRACK157 
TRACK158 
TRACK159 
TRACK160 
TRACK161 
TRACK162 
TRACK163 


ci  rj  o  o  o  o  o 


PRINT  *,  'TO  HAVE  CODE  COMPUTE  CONE  INDEX.  INPUT  CI-0.' 

PRINT  *.  'TO  HAVE  CODE  COMPUTE  ROLLING  RESISTANCE.  INPUT  SF*0. 
PRINT  *.  'INPUT  CI . A.SM.SN.SXI .CD.CLAMDA.PHI .G.SF.ETA.CHI' 

READ  *.  CI. A.SM.SN.SXI. CD. CLAMDA.PHI .G.SF. ETA, CHI 

C  =  A  -  SM 

IF  (SXI  .LT.  1.)  C  =  0. 

IF  (CI  .LE.  0.)  CALL  CONE 

IF  (SF  .LE.  0.)  SF  =  .045  +  (2.3075/(CI-VCIl+6.5) ) 

PRINT  *, 

PRINT  *,  'INPUT  RUNTYPE  :  ' 

PRINT  *,  '  FOR  STEADY-STATE  RUN.  INPUT  0' 

PRINT  *,  '  FOR  TRANSIENT  RUN,  INPUT  1' 

PRINT  *,  'INPUT  YOUR  CHOICE' 

READ  *,  IRUN 

IF  (IRUN  .EQ.  0)  CALL  STEADY 
IF  (IRUN  .EQ.  1)  CALL  TRANS 

SYSTEM  DEPENDENT  ROUTINE  TO  GET  PROCESSOR  TIME 

CALL  PTIME(PTU) 

PTU  =  (PTU-PTI )*60 
PRINT  250,  PTU 
PRINT  246 


FORMAT  < A40) 

FORMAT ( IH1 ) 

FORMATS, '  PROCESSOR  TIME  (CPU-MIN) 


=  ' » 1PE13 ,5) 


************************************************************** 

*  * 

*  THIS  SUBROUTINE  IS  FOR  STEADY-STATE  TURNING  MOTION.  * 

*  * 
************************************************************** 

SUBROUTINE  STEADY 

PARAMETER  (  N2  =  100, N4  =  40  ) 

REAL  L 

LOGICAL  LOGRE.LSTB.LVX 

CHARACTER  TESTN*40,DATE*8 

CHARACTER  LAB2*5, INPUT, VFMT*3»LFMT*64< 4) 

DATA  VFMT  /3H(V)/ 

COMMON  /ARRAYS  /  TP(N2) »VXP(N2) , VYP(N2) ,FCXP(N2> »FCYP(N2) » 
t  VSX1P(N2) ,VSX2P(N2) »R0P(N2) ,PTEP(N2>  »PTSP(N2) , 

t  DXTP<N2)»DYTP(N2),MDP(N2),WP(N2)»VX1P(N2)» 

1  VX2P(N2) ,EP(N2) ,VP(N2) ,PP(N2) »PT1P(N2) »PT2P(N2) 


TRACK164 
TRACK165 
TRACK166 
TRACK167 
TRACK168 
TRACK169 
TRACK170 
TRACK171 
TRACK172 
TRACK! 73 
TRACK174 
TRACK175 
TRACK176 
TRACK177 
TRACK178 
TRACK179 
TRACK180 
TRACK181 
TRACK182 
TRACK183 
TRACK184 
TRACK185 
TRACK186 
TRACK187 
TRACK188 
TRACK189 
TRACK190 
TRACK191 
TRACK192 
TRACK193 
TRACK 194 
TRACK195 
TRACK196 
TRACK197 
TRACK198 
STEAD001 
STEAD002 
STEAD003 
STEAD004 
STEAD005 
STEAD006 
STEAD007 
STEAD008 
STEAD009 
STEAD010 
STEAD011 
STEAD012 
STEAD013 
STEAD014 
STEAD015 
STEAD016 
STEAD017 
STEAD018 
STEAD019 
STEAD020 
STEAD021 


I 

$ 

kv 


I 


KK 


$ 


COMMON  /INOUT  /  E.VX1.VX2.V.RIOL.LOGRE.LVX 


COMMON  /IN  / 
COMMON  /VPARAM/ 


DTIM1.DTIM2.  UPRV . WDPRV » VXPRV. VYPRV  » TESTN.DATE 
WT  .L.H.D.B.SL.CX.THETAA.THETAD. IZ.VCI1 


I 


COMMON  /3PARAM/  Cl .  A.SM. SN. C.SXI .CD.CLAMDA.PHI » G.SF » ETA.CHI 
COMMON  /OUT  /  FCX  »FCY » VSX1 » VSX2. PTE. PTS. WD. W» VX» VY . 
l  P.Z1.Z2.IDIR.VXD.VYD.PT1.PT2 

COMMON  /MISC  /  GRAV.DRAD. VMP.RDEG 

PRINT  *» 

PRINT  *.  'INPUT  CASE*.  1.2.3. 4.  OR  ?  ' 

READ  *.  INPUT 

IF  (INPUT  .EQ.  1H?)  GO  TO  900 
DECODE< INPUT»VFMT)  ICASE 
IF  (ICASE. GE.l  .AND.  ICASE. LE. 4)  GO  TO  915 
900  PRINT  910 

PRINT  245.  LFMT 
PRINT  916 
READ  *.  ICASE 

910  FORMAT ( '  CASE  SELECTIONS! ' ) 

DATA  LFMT  / 

i  64H  CASE  C13 . VARY  V  (VELOCITY) 

t  INPUT  E  (STEERING  RATIO)  , 

S  64H  CASE  C2] . VARY  V  (VELOCITY) 

1  INPUT  RO  (TURNING  RADIUS). 

%  64H  CASE  C33 . . . .VARY  E  (STEERING  RATIO) 

t  FIXED  V  ( VELOCITY )  . 

%  64H  CASE  C43 . . .VARY  RO  (TURNING  RADIUS) 

%  FIXED  V  (VELOCITY)  / 


916  FORMAT ( '  INPUT  YOUR  CHOICE! 
915  PRINT  *» 

IF  (ICASE  .GT.  2)  GO  TO  5 
IF  (ICASE  .EQ.  2)  GO  TO  4 

PRINT  *.  'INPUT  E.DV.HPAV  ' 
READ  *.  E.DV.HPAV 
LOGRE  =  .FALSE. 

GO  TO  75 

4  CONTINUE 
LOGRE  =  .TRUE. 

PRINT  *,  'INPUT  RO.DV.HPAV 
READ  *.  RO.DV.HPAV 
RI  =  RO 

RIOL  =  12.4RI/L 

75  CONTINUE 
V  =  0. 

NP  =  N4 
IHP  =  0 

IF  (HPAV  .LE.  0.)  IHP  =  1 
GO  TO  7 

5  CONTINUE 
DV  =  0. 

IHP  =  1 


STEAD022 

STEAD023 

STEAD024 

STEAD025 

STEAD026 

STEAD027 

STEAD028 

STEAD029 

STEAD030 

STEAD031 

STEAD032 

STEAD033 

STEAD034 

STEAD035 

STEAD036 

STEAD037 

STEAD038 

STEAD039 

STEAD040 

STEAD041 

.STEAD042 

STEAD043 

.STEAD044 

STEAD045 

.STEAD046 

STEAD047 

JSTEAD048 

STEAD049 

STEAD050 

STEAD051 

STEAD052 

STEAD053 

STEAD054 

STEAD055 

STEAD056 

STEAD057 

STEAD058 

STEAD059 

STEAD060 

STEAD061 

STEAD062 

STEAD063 

STEAD064 

STEAD065 

STEAD066 

STEAD067 

STEAD068 

STEAD069 

STEAD070 

STEAD071 

STEAD072 

STEAD073 

STEAD074 

STEAD075 

STEAD076 

STEAD077 


75 


o  n  o  o  o 


IF  (ICASE  .EQ.  4)  GO  TO  6 
PRINT  165 »  N2 

PRINT  *.  'INPUT  V.E.DE.NP' 

READ  *»  V.E.DE.NP 
LOGRE  =  .FALSE. 

DR  =  0. 

GO  TO  7 

6  CONTINUE 
PRINT  165.  N4 

PRINT  *.  'INPUT  V.RO.DR.NP' 

READ  *»  V.RO.DR.NP 
LOGRE  *  .TRUE. 

DE  =  0. 

RI  =  RO 

RIOL  =  12.IRI/L 

7  CONTINUE 
PRINT  246 
PRINT  250 

PRINT  245.  LFNT < ICASE ) 

PRINT  207 

IF  (ICASE. EQ.l  .OR.  ICASE. EQ. 2) 
l  PRINT  202.  DV.HPAV.TESTN.DATE 

IF  (ICASE  .EQ.  3)  PRINT  203.  DE.NP.TESTN.DATE 
IF  (ICASE  .EQ.  4)  PRINT  204.  DR.NP.TESTN.DATE 
PRINT  170 

PRINT  180.  WT.L.H.D.B.SL.CX.THETAD.THETAA.IZ.VCIl 
PRINT  *. 

PRINT  190 

PRINT  200.  CI»A»SM»C»SN.IFIX(SXI>.CD» 

S  CLAMPA.PHI.G.SF.ETA.CHI 

PRINT  207 
PRINT  205 
PRINT  209 

LSTB  =  .FALSE. 

V  =  0*17.6 
DV  =  DV*17.6 
DR  =  12.*DR/L 
DTIN1  =  0. 

WDPRV  =  0. 

VXPRV  =  0. 

VYFRV  =  0. 

LVX  =  .FALSE. 


STEADY-STATE  CALCULATIONS 


DO  400  II  =  l.NP 
11  CALL  AGIL(*115) 

IF ( 1 1 . EQ .  1  .OR.  ICASE. GT. 2)  GO  TO  390 
IF < <P-PP( II-l ) )  .GT .  0.1)  GO  TO  115 


STEAD078 
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m 


E; 


ft 


m 


390  CONTINUE 
PP<II)  =  P 
ZIP  =  21 
Z2P  =  Z2 
VP(II)  =  V4UHP 
ROP< 1 1 )  =  RI0L*L/12. 

EP  < II)  =  E 

IF  <E  .LT.  1.)  EP< II )  =  -l./E 
V  =  V+DV 

VSXIP(II)  =  VSX1IVMP 
VSX2P (II)  =  VSX2*VMP 
VX2P< I I )  =  VX2CVMP 

vxipui)  =  vxuvHP 

WDP(II)  =  WD*RDEG 


PTEP (II)  =  PTE 
PTSP(II)  =  PTS 
PTIP(II)  =  PT1 
PT2P( II )  =  PT2 
FCXP(II)  =  FCX 
FCYP<II)  =  FCY 

PRINT  14.  1 1  »F'P<  II )  »EP(  1 1 ) »ROP( II ) »VP(II)»VX1P(II) » VX2P( I I)r 
t  VSX1P<II>,  VSX2P( 1 1 ) « WBP<  I I)iFCXP(II) , FCYP( 1 1 ) » 

S  PTEP(II),PTSF<II) 

392  CONTINUE 

IF  < ICASE  .GT.  2)  GO  TO  160 
IF<P.GT. 0.555  .AND.  LSTB)  GO  TO  404 

CALCULATES  MAXIMUM  POWER 

IF(IHP)  393.393,396 

393  IF( PTEP< 1 1 )-HPAV)  396,394,394 

394  PFC  =  (PTEP(II)-HPAV)  /  <  PTEP<  1 1  )-F'TEP<  1 1-1 ) ) 

IHP  =  1 

ROPF  =  ROP(II)  -  (ROP(II)-ROP(II-1))*PFC 
YPF  =  UP(II)  -  (VP(II)-VPUI-l)  )*PFC 
VX1F  =  VX1P<II)  -  <VX1P<II)-VX1P<II-1))*PFC 
VX2F  =  VX2P( II )  -  (VX2F(II)-VX2P(II-1))*PFC 
VSX1F  =  VSXIP(II)  -  (MSX1P(  II  )-VSXlP(  1 1  —  1  >  )#F’FC 
VSX2F  =  VSX2PC II)  -  ( VSX2F( I I )  -VSX2P< 1 1- 1 ) )*PFC 
PTEF  =  PTEP(II)  -  ( PTEP  <  1 1 )  -F‘TEP(  1 1  —  1 ) )  <PFC 
PTSF  =  PTSP< II )  -  (PTSP(II)-PTSP(II-l) >*PFC 
FCXF  =  FCXF(II)  -  (FCXP(II)-FCXP(II-l)) *PFC 
FCYF  =  FCYF'UI)  -  <FCYF<  II  >-FCYP(Il-l )  )*PFC 
PF  =  PP(II)  -  < PP til) -PP< 1 1  —  1 ) )*PFC 
EF  =  EP(II)  -  (EF'( II  >-EP( 1 1-1 )  )*PFC 
WDF  =  UDP(II)  -  < MDP< II) -WDP< 1 1-1 ) ) *PFC 
PRINT  *,  'POWER  LIMIT' 

PRINT  14,  1 1 »PF » EF »ROPF , VPF , VX1F, VX2F , VSX1F , VSX2F , 


WDF, FCXF, FCYF, PTEF, PTSF 


PRINT  *, 


396  CONTINUE 

IF(LSTB)  GO  TO  400 


STEAD134 
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STEAD165 
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STEAD168 
STEAD169 
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STEAD175 
STEAD176 
STEAD177 
STEAD178 
STEAD179 
STEAD180 
STEAD181 
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STEAD183 
STEAD184 
STEAD185 
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\ 

!  . 
* 

A 

397 

398 

160 

400 

IF(II  .£0.  1)  GO  TO  400 

IF<P“.5>  397.462.460 

IF < ( VSX2P (II) -VSX2P < 1 1 - 1 ) ) /DV/VNP-1 . )  398 » 398 » 464 

1F< (VSX1P<I1 >-VSXlP< 11-1) )/DV/VMP-l .  >  400.400.466 

CONTINUE 

E  =  E+DE 

RIOL  »  RIOL  -  DR 

CONTINUE 

STEAD190 

STEAD191 

STEAD192 

STEAD193 

STEAD194 

STEAD195 

STEAD196 

STEAD197 

STEAD198 

STEAD199 

Jj  c 

altnUiVv 

{j  c 

STEAD201 

s« 

IF  (ICASE.GT.2)  GO  TO  480 

STEAD202 

Jj 

PRINT  165*  N4 

STEAD203 

? 

IF  < .NOT >  LSTB)  GO  TO  480 

STEAD204 

1  c 

STEAD205 

1 

404 

PRINT  207 

STEAD206 

2 

PRINT  *.  •  STABILITY  LIMIT  0NSLAB2 

STEAD207 

5 

PRINT  t, 

STEAD208 

i 

PRINT  14.  N3 . PS » ES » ROPS » VPS » VX 1 S » VX2S » VSX1 S . 

STEAD209 

1  VSX2S.WDS.FCXS.FCYS.PTES.PTSS 

STEAD210 

GO  TO  99 

STEAD211 

1  c 

STEAD212 

460 

N3  =  II 

STEAD213 

PDC  =  (PP(N3)-.5)  /  <PP(N3)-PP(N3-1>) 

STEAD214 

> 

IF<P  .GT.  0.555)  LSTB  =  .TRUE. 

STEAD215 

• 

GO  TO  463 

STEAD216 

'•  c 

STEAD217 

i 

462 

N3  =  II 

STEAD218 

1 

PDC  *  0. 

STEAD219 

»* 

463 

LAB2  =  '  P  ' 

STEAD220 

GO  TO  468 

STEAD221 

t  c 

STEAD222 

tt 

464 

N3  =  II-l 

STEAD223 

) 

PDC  =  0. 

STEAD224 

1 

LAB2  =  '  VS2' 

STEAD225 

1 

GO  TO  468 

STEAD226 

c 

STEAD227 

466 

N3  =  II-l 

STEAD228 

• 

PDC  =  0. 

STEAD229 

* 

LAB2  =  '  VS1' 

STEAD230 

i 

GO  TO  468 

STEAD231 

1  C 

STEAD232 

1 

» 

467 

N3  =  II  -  1 

STEAD233 

: 

LSTB  =  .TRUE. 

STEAD234 

: 

PDC  =  0. 

STEAD235 

LAB2  =  '  DV  ' 

STEAD236 

c 

STEAD237 

468 

CONTINUE 

STEAD238 

ROPS  =  R0P<N3)  -  «R0P(N3)-R0P<N3-1))*PDC 

STEAD239 

VPS  =  VP<N3)  -  <VP<N3)-VP(N3-1>)*PDC 

STEAD240 

VSX1S  =  VSX1P<N3)  -  <VSX1P(N3)-VSX1P(N3-1))*PDC 

STEAD241 

k 

VSX2S  =  VSX2P(N3)  -  (VSX2P(N3)-VSX2P<N3-l) )*PDC 

STEAD242 

t 

VX1S  =  VX1P<N3>  -  <VX1P<N3>-VX1P(N3-1))*PDC 

STEAD243 

VX2S  =  VX2P(N3)  -  <VX2P(N3)-VX2P(N3-1 ) )*PDC 

STEAD244 

PTES  =  PTEP(N3)  -  (PTEP(N3)-PTEP(N3-1))*PDC 

STEAD245 

78 

■  ..... 

PTSS  =  PTSP<N3)  -  <PTSP(N3)-PTSP(N3-1))*PDC 

FCXS  =  FCXP<N3)  -  <FCXP<N3)-FCXP<N3-1>>*PBC 

FCYS  =  FCYP<N3)  -  (FCYP<N3)-FCYP(N3-1 ) )*PDC 

PS  =  PP(N3)  -  (PP(N3)-PP(N3-1))*PDC 
ES  =  EP(N3)  -  <EP(N3)-EP(N3-1))*PDC 
UDS  =  WBP(N3)  -  <UDP(N3)-UDP<N3-1))*PBC 
IF(LSTB)  GO  TO  404 
LSTB  =  .TRUE. 

GO  TO  400 

115  IFUI.EQ.l  .OR.  ICASE.GT.2)  GO  TO  11? 

IF(DV  .LT.  1.)  GO  TO  467 
DV  =  DV/2. 

V  =  V  -  DV 
P  =  PP<II-1) 

Z1  =  ZIP 
Z2  =  Z2P 
GO  TO  11 


'THE  VEHICLE  HILL  NOT  GO  IN  THIS  PARTICULAR  SOIL.' 


119  PRINT  *,  ' 
GO  TO  99 

480  PRINT  207 

99  RETURN 


14  FORMAT (13X,I3»1X,F6.2»F7.2,F10.2,1X,F5.1»2X,F5.1,1X»F5.1»2F7.2» 
t  F7 . 1 »2F7 . 2»2F7 . 1 ) 

165  FORMAT ( '  MAXIMUM  POINTS  NP  SET  TO'. 15) 

170  FORMAT <2X. 'V)T(IN)'»3X»'L(IN)'f3X»'H(IN)'.2X.'B(IN)'»4X» 

*  'B<IN)'.3X»'SL<IN)'.1X.,CX(IN)'.4X. 'THETAD(DEG) ' .3X» 

l  ' THETAA(DEG) ' . IX. ' IZ< IN  LB/SECM2) '  ,3X, 'VCIl' > 

180  FORMAT (1X,F8.1»2X»F5.1,2(3X»F4.1),4X»F5.1»5X,F4.1,2X»F5.1, 

(  7X.F5.1»10X.F5.1»7X»F8.1f8X.F4.1) 

190  FORMAT <4X» ' Cl' .5X» ' A(PSI ) ' . IX. 'SM<PSI ) ' . IX . ' C (PSI ) ' . IX, 

S  '  SN( 1/PSI ) ' ,2X , ' SXI ' , 3X, 'CD<  PSI ) ' , IX 

I  '  CLAMDA< SEC/IN) '»4X» 'PHI (BEG)', 2X, 'G(PSI/IN> ' ,2X» 

l  'SF' ,4X, ' ETA < BEG) ' ,1X, 'CHI (DEG) ' ) 

200  FORMAT (2X»F5.1»1X»2(3X»F5.1)»3X»F4.1»3X»F6.4,5X,I1,5X»F5.1, 
i  6X,F6.3,10X,F4.1,6X,F5.1,3X,F5.3»4X,F4.1,4X,F5.1) 

202  FORMAT ( 2X» ' DV  =  ' » F6 . 1 » ' ( MPH ) ' , 1 X » ' HPAV  =  ' ,F6. 1 , ' (HP) ' » 

S  3X , ' TEST  NO.  =  '  .A40.3X.A8.3X.//) 

203  FORMAT < 2 X , 'DE  =  ' .F6.2.1X, 'NP  =  ',13, 

S  6X , ' TEST  NO.  =  ' ,  A40.3X.A8.3X,//) 

204  FORMAT < 2X, 'DR  =  ' ,F6. 1 » ' (FT) ' , IX, 'NP  =  ',13, 

l  6X, ' TEST  NO.  =  ' ,A40,3X, A8.3X,//) 

205  FORMAT (14X,'PT',4X,'P/L',4X,'E',6X,'R0',7X, 'VEL' ,4X, 'VX1 ' , 

t  3X,'VX2',3X,'VS1',4X,'VS2',4X,'WD',5X,'FCX', 

t  4X , ' FCY ' , 4X » 'PTE ' ,  4X, ' F’TS  ' ) 

207  FORMAT </,lX,119(lH=) ,//) 

20?  FORMAT (34X, 'FT ' , 7X , ' MPH ' , 4X , 2 ( ' MPH ' ,3X ) , ' MPH ' , 4X, ' MPH ' ,2X , 

*  '  DEG/SEC ',3X»'G',6X»'G',1X»2<5X,' HP '),//) 

245  FORMAT (28X.A64) 

246  FORMAT(lHl) 

250  F0RMAT<46X, 'STEADY-STATE  TURNING  MOTION') 
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END 

************************************************************** 
*  * 

*  THIS  SUBROUTINE  IS  FOR  TRANSIENT  MOTION.  * 

*  * 
************************************************************** 

SUBROUTINE  TRANS 

PARAMETER  <  N2  =  100  ) 

REAL  L> IZ 

LOGICAL  LOGRE,LVX,LSTB,LOGV 


CHARACTER  TESTN*40 , DATE*8 
CHARACTER  INPUT, VFMT*3, LFMT*76(4> 
DATA  VFMT  /3H(V>/ 


COMMON  /ARRAYS/ 
l 
i 
l 

COMMON  /INOUT  / 
COMMON  /IN  / 
COMMON  /VF'ARAM/ 
COMMON  /SPARAM/ 
COMMON  /OUT  / 
l 

COMMON  /MISC  / 


TP<  N2) ,  VXF'<  N2) » VYP(N2> , FCXP (N2) ,FCYP (N2) , 
VSX1F(N2) ,VSX2P(N2),R0P(N2),PTEP(N2) ,PTSP(N2>, 
DXTP(N2) ,DYTP(N2) ,WDP(N2) ,WP<N2) ,0X1P(N2) , 
VX2P(N2) »EP(N2) ,VP<N2) »PP(N2) »PT1P(N2) »PT2P(N2) 
E,VX1,VX2,V,RI0L, L06RE ,LVX 

DTIM1 , DTIM2, WPRVfWIiPRV, VXPRV, VYPRV, TESTN, DATE 
UT,L,H,D,B,SL,CX,THETAA,THETAD,I2,UCI1 
CI,A,SM,SN,C,SXI,CD,CLAMDA,PHI,G,SF,ETA,CHI 
FCX,FCY, VSX1 , VSX2,PTE,PTS,WD,W,VX,VY, 

F',Z1  ,Z2,  IDIR,VXD»VYD,PT1  »PT2 
GRAVfDRAD,  VMF'jRDEG 


THE  DIGITIZED  TIME  HISTORY  (FOR  CASE  3)  AND  THE  FUNCTION 
(FOR  CASES  1,2,  AND  4)  USED  FOR  SAMPLE  PROBLEMS 

DATA  EP  /10*1. 1,10*1. 0»10*-1. 1,70*1./ 

EXPS(X)  =  EXF'( AMAX1  (X,-88. ) ) 

FT( T>  XI » X2»  X3’X4 )  = 

l  XI  +  X2* ( EXPS ( -X3*T )  +  X3*T*EXPS(-X4*T)  -  1.) 


note:  to  drive  the  MODEL,  THE  above  may  be  REPLACED  BY  ANY 
DIGITIZED  TIME  HISTORY  AND/OR  ANY  FUNCTION 


PRINT  *, 

PRINT  165,  N2 

PRINT  *,  'INPUT  DT,NP, IPRINT ' 

READ  *,  DT,NP, IPRINT 

PRINT  *, 

PRINT  *,  'INPUT  CASE:  1,2, 3, 4,  OR  ?  ' 
READ  *,  INPUT 

IF  (INPUT  .EQ.  1H?)  GO  TO  900 
IF  (INPUT  .EQ.  1H1 )  ICASE  *-  1 


STEAD302 
STEAD303 
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TRANS015 
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TRANS053 
TRANS054 
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c 


c 


IF  (INPUT  .EQ.  1H2 )  ICASE  =  2  TRANS055 

IF  (INPUT  .EQ.  1H3)  ICASE  =  3  TRANS056 


IF  (INPUT  .EQ.  1H4 )  ICASE  =  4  TRANS057 

IF  (ICASE. GE.l  .AND.  ICASE. LE. 4)  GO  TO  915  TRANS058 

900  PRINT  910  TRANS059 

PRINT  245*  LFMT  TRANS060 

PRINT  916  TRANS061 

READ  t>  ICASE  TRANS062 

910  FORMAT ( / '  CASE  SELECTIONS! ' )  TRANS063 

DATA  LFMT  /  TRANS064 

S  54H  CASE  m . INPUT  V.E  (VELOCITY  *  STEERING  RATIO)  *  TRANS065 

1  54H  CASE  C2I . INPUT  V»RI  < VELOCITY  *  TURNING  RADIUS)  *  TRANS066 

S  76H  CASE  C3I... INPUT  E.VI  (STEERING  RATIO* INITIAL  VELOCITY)  *TRANS067 
1  VARY  V  (VELOCITY)  *  TRANS068 

S  54H  CASE  14] . INPUT  VX1.VX2  (TRACK  VELOCITIES)  /  TRANS069 

916  FORMAT (  'INPUT  YOUR  CHOICE!  ')  TRANS070 

915  PRINT  *»  TRANS071 

TRANS072 


LOGRE  =  .FALSE. 

LVX  =  .FALSE. 

IF  (ICASE  .GT.  2)  GO  TO  3 

PRINT  t,  'INPUT  HPAV*V1*V2*V3*V4' 

READ  **  HPAV* VI * V2* V3* V4 
VI  =  VI 

IF  (ICASE  .EQ.  2)  GO  TO  2 
PRINT  *,  'INPUT  El  *  E2  *E3  *E4  ' 

READ  *»  El »E2f E3*E4 
GO  TO  7 

NOTE!  CASE  2  COULD  BE  MODIFIED  TO  INPUT  TRAJECTORY 

INSTEAD  OF  RI 

2  PRINT  *»  'INPUT  RI1.RI2*RI3»RI4' 

READ  RI1 «RI2»RI3»RI4 
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LOGRE  =  .TRUE. 

GO  TO  7 

3  IF  (ICASE  .EQ  .4)  GO  TO  4 
PRINT  **  'INPUT  HPAViVI' 

READ  *»  HPAVf VI 

GO  TO  7 

4  LVX  =  .TRUE. 

PRINT  **  'INPUT  VXllfVX12»VX13fVX14' 
READ  **  VXllrVX12fVX13fVX14 
FRINT  *i  'INPUT  VX21»VX22*VX23»VX24' 
READ  *f  VX21fVX22fVX23»VX24 
VI  =  AMAX1 ( VX1 1  *  VX21 ) 

7  PRINT  246 
PRINT  250 

PRINT  245*  LFMT ( ICASE ) 

IF  (ICASE  .EQ.  2)  PRINT  225 
PRINT  207 

PRINT  203*  BTiNP.TESTN.BATE 
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IF  (ICASE  .EQ.  1)  PRINT  210*  HPAV. VI . V2* V3. V4 *E1 .E2.E3.  E4 

IF  (ICASE  .EQ.  2)  PRINT  220*  HPAV*V1 »V2»V3* V4.RI1 .RI2.RI3.RI4 

IF  (ICASE  .EQ.  3)  PRINT  230*  HPAV.VI 

IF  (ICASE  .EQ.  4)  PRINT  240*  VI  * VX1 1  * VX12* VX13* VX14* 

%  VX21* VX22.VX23. VX24 

PRINT  ** 

PRINT  170 

PRINT  180*  UT*L  »H*D*B  *SL  *CX  *  THETAD*THETAA* IZ  * VCI 1 
PRINT  ** 

PRINT  190 

PRINT  200*  CI.A*SN*C.SN.IFIX(SXI).CD. 

%  CLAHDA.PHI*G*SF,ETA»CHI 

PRINT  207 
FRINT  205 
PRINT  209 
C 

C  SET  UP  INITIAL  CONDITIONS 
C 

T  =  DT 

VX  =  VI*17.6 
VY  =  0. 

UD  =  0. 

DXTPS  =  0. 

DYTPS  =  0. 

THETA  =  0. 

L0GV  =  .FALSE. 

DTIM1  =  1 ./DT 
DTIH2  =  DT/2. 

DT012  =  DTIH2/12. 

DXTPRV  =  VX 
DYTPRV  =  VY 
THDF'RV  =  0. 

UPRV  =  0. 

UDPRV  =  UD 
VXPRV  =  VX 
VYPRV  =  VY 


TRANSIENT  CALCULATIONS 


DO  400  II  =  l.NP 
C 

LSTB  =  .FALSE. 

DA  =  8.8 
DV  =  8.8 

GO  TO  (10*20*30*40)  .ICASE 
C 

10  CONTINUE 

E  =  FT(T.E1*E2.E3*E4) 

V  =  FT(T*V1*V2.V3>V4)*17,6 

IF(LOGV)  V  =  AHIN1(VP(II-1)*17.6  +  2.*DA*DT  *  V) 
GO  TO  11 
C 

20  CONTINUE 

RIOL  =  FT(T*RI1.RI2*RI3*RI4)*12./L 
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V  =  FT(T.V1»V2.V3.V4)*17.6 

IF<L06V>  V  =  AMIN1(VP<II-1>*17.6  ♦  2.*DA*DT  »  V) 
60  TO  11 

CONTINUE 

V  =  SQRT<VXPRV**2+VYPRV**2>  +  DA*DT 
E  =  EP(II) 

IF ( .NOT. LOGV)  GO  TO  11 
El  =  E 

IF( VPREV  .  GE.  VI)  GO  TO  48 
IF<ABS<EP< 1 1-1 )-l . )  .LE.  0.05)  GO  TO  47 
IF( VPREV  .LT.  VI)  E  =  EP(II-l)  -  .05BIDIR 
GO  TO  49 

E  =  1.  *  <EI-1.)/ABS<EI-1. >*.000001 
CONTINUE 

IF (RIOLME/EI-l • ) )  11.48*46 
E  =  El 

LOGV  =  .FALSE. 

GO  TO  11 

CONTINUE 

VX1  =  FT<T.VXll.VX12»VX13tVX14)*17.6 
VX2  =  FT ( T » VX21 » VX22i VX23» VX24)#17»6 
E  =  VX1/VX2 

CALL  AGIL0115) 

ZIP  =  Z1 
Z2P  =  Z2 
PP< II )  =  P 

IF(LVX)  V  *  SQRT<VX**2  +  VY**2> 

VP<II)  =  V*VMP 
VPREV  =  V 

TF (LVX  .OR.  HPAV.LE.O.)  GO  TO  160 
RATIO  =  ABS<PTE/HPAV) 

IF<RATIO  .GT.  1.0)  GO  TO  154 
IF (RATIO  . GE.  0.95)  GO  TO  160 
IF<ICASE.NE.3>  GO  TO  156 
IF<LSTB)  GO  TO  160 

V  =  V+DAIDT 
GO  TO  11 

CONTINUE 

V  =  V-DA4DT 

DA  =  AMAX1 ( DA/1 . 2. . 05 ) 

IF(V.LE,VI  .AND.  ICASE.EQ.3)  LOGV  =  .TRUE. 

IF ( ICASE  .LT.  3)  LOGV  =  .TRUE, 

GO  TO  11 

LOGV  =  .FALSE. 

CONTINUE 

UDP(II)  =  UD4RDEG 

ALPHD  =  <VX*VYD-VY*VXD)  /  V**2 

THED  =  WD  -  ALPHD 

ROP(II)  =  V/THED/12. 
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THETA  =  THETA  +  ( THED+THDPRV ) *DTIM2 
DXT  =  V*C0S ( THETA) 

DYT  =  V*SIN(THETA> 

DXTPS  =  DXTPS  +  ( DXTPRV+DXT ) IDT012 

DYTPS  =  DYTPS  t  <DYTPRV+DYT>*DT012 

DXTP< II )  =  DXTPS 

DYTP(II)  =  DYTPS 

DXTFRV  =  DXT 

DYTF'RV  =  DYT 

THDPRV  =  THED 

WDPRV  =  WD 

VXPRV  =  VX 

VYFRV  =  VY 

HPRV  =  U 

UP(II)  =  U*RDEG 

Y  <  1 1 )  =  T 

E  P  < II)  =  E 

IF  (E  .LT.  1.)  EP(II)  =  -l./E 

FCXP(II)  =  FCX 

FCYP(II)  =  FCY 

VXIP(II)  =  VX1JVHP 

VX2P< II )  =  VX2*VMP 

VSXIP(II)  =  VSX14VMP 

VSX2P(II)  =  VSX2*VMP 

VXP(II)  =  VX*VMF 

VYP(II)  =  VY*VHP 

PTEP(II)  =  PTE 

FTSP(II)  =  PTS 

P  T 1 P  < 1 1 )  =  PT1 

PT2P( II )  =  PT2 

T  =  T  +  DT 

IF  <  HOD( II , IPRINT) ,NE .0 )  GO  TO  400 
PRINT  14,  1 1 ,  TF  ( 1 1 )  >  F'P  ( 1 1 ) » 

l  EP( I I ) , ROP( II ) , VP(II)»VX1P<II),VX2P(II), 

*  VSX1P( II ) • VSX2P< II ) » UPC  II ) , WDP( II ) ,  FCXP(II), 

i  FCYP(II),  FTEP  ( 1 1 )  ,PTSP<  II ) ,  DXTF'(  II ) ,  DYTP(  1 1 ) 

C 

400  CONTINUE 

C  - - - - 

c 

PRINT  207 
GO  TO  9? 


115  IF 

(II 

.EQ.  1) 

GO 

TO 

119 

0  = 

0 

-  DV 

IF  ( 

DV 

.LT.  1.) 

GO 

TO 

116 

DV  =  DV/2. 

LSTB  =  .TRUE. 

P  =  PP<II-1) 

Z1  =  ZIP 
Z2  =  Z2P 

IF(.NOT.LVX)  GO  TO  11 
0X1  =  VX1P( I 1-1 )*17 . 6 
VX2  =  VX2P ( II-l )*17.6 
E  =  VX1/VX2 
GO  TO  11 
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116  PRINT  207 

PRINT  *»  'VEHICLE  GOES  UNSTABLE  AT  THIS  POINTS 
t  '  (TOO  VIOLENT  MANEUVER  ATTEMPTED)' 

PRINT  *» 

II  =  II-l 

PRINT  1 4 »  II »TP( II ) »PP( II) » 

t  EP  < 1 1 ) *ROP( II)»VP(II).VX1P(II).VX2P(II). 

4  VSX1P( II ) »VSX2P< I X ) » UP (II) .UDPdl ) »  FCXP(II). 

4  FCYP(II),  PTEP( II ) »PTSP< II ) »DXTP( II) . DYTP( II ) 

GO  TO  99 
C 

119  PRINT  *.  'THE  VEHICLE  HILL  NOT  GO  IN  THIS  PARTICULAR  SOIL.' 

C 

99  RETURN 

C  ============================================================== 

c 

14  FORMAT <1X»I3»F6.2»F7.2»F6.2*F10.2»3F6.1*2F7.2»F7.1*F6. 1» 

*  2F6 .2 >2F7 . 1 i2F8 . 1 ) 

165  FORMAT ( '  MAXIMUM  POINTS  NP  SET  TO'. 15) 

170  FORMAT (2X i ' WT( IN) ' » 3X» 'L< IN) ' .3X. 'H< IN) ' »2X» 'D( IN) ' » 4X. 

&  'B(IN)'»3X»'SL(IN)'.1X»'CX(IN)'»4X> 'THETAD(DEG) ' »3X» 

4  'TH£TAA(DEG)'»1X»'IZ(IN  LB/SEC<*2 ) ' » 3X » ' VC  1 1 ' ) 

180  F0RHAT(lX»F8.1»2X»F5.1»2<3X»F4.t>»4X.F5.1.5X»F4.1»2X.F5.1» 

4  7X»F5.1»10X>F5.1»7X»F8.1»8X.F4.1) 

190  F0RHAT(4X.'CI'»5X.'A(PSI>'»1X. 'SM(PSI > ' » IX. 'C(PSI) ' . IX. 

4  'SN< 1/PSI ) ' »2X. 'SXI' .3X» 'CD(PSI) ' »1X 

4  'CLAMDA( SEC/IN) '  .4X» 'PHI (DEG) ,.2X.'G(PSI/IN),f2X. 

4  'SF',4X.  'ETA(DEG) '  »1X.  'CHI (DEG) ' ) 

200  F0RNAT(2X.F5.1.1X.2<3X,F5.1)»3X»F4.1»3X»F6.4.5X»I1»5X,F5.1. 

4  6X»F6.3>10X.F4.1»6X»F5.1.3X»F5.3»4X>F4.1»4X»F5.1> 

203  FORMAT ( 2X. ' DT  =  ' »F4.1 . ' (SEC) ' »3X. 'NP  =  '.13. 

4  12X. 'TEST  NO.  =  ' .A40.3X.A8.3X.//) 

205  FORMAT (2X.'PT',3X.'T'.5X.'P/L'.4X.'E'.6X.'R0'.7X.'VEL'.3X,'VX1'. 

*  3X»'VX2'.3X»'VS1'.4X. 'VS2'.5X.'U'.5X.'WD'.4X.'FCX'. 

4  3X.'FCY',3X»'PTE'.4X. 'PTS'.5X.'DXT'.5X.'DYT'> 

207  FORMAT </. IX. 119 ( 1H=) »//) 

20?  FORMAT ( 6X » 'SEC '»18X»'FT'»7X»3( 'MPH' . 3X ) » 2( 'MPH' .4X) . 'DEG' . 

4  2X. ' DEG/SEC '.2X» 'G' .5X. 'G' .2(5X» 'HP' ) .2(6X» 'FT' ) .//) 

210  FORMAT ( 2X . 'HPAV  =  ' .F6 . 1 ,' (HP) ' »7X» ' DRIVER  COEFFICIENTS  FOR'. 

4  '  VELOCITY  5'»F6.2»3('»'.F6.2)» /»50X» 

4  'STEERING  RATIO* '»F6.2.3('.'.F6.2)> 

220  F0RMAT(2X. 'HPAV  =  ' »F6.1. ' (HP) ' »7X» 'DRIVER  COEFFICIENTS  FOR'. 

4  '  VELOCITY  J '»F8.1.3< '.' .F8.1)./»50X» 

4  'TURNING  RADIUS* '»F8.1»3('.'»F8.1>) 

225  FORMAT(25X.'NOTEJ  THIS  CASE  COULD  BE  MODIFIED  TO  INPUT'. 

4  '  TRAJECTORY  INSTEAD  OF  RI') 

230  F0RMAT(2X. 'HPAV  =  ' ,F6. 1 »4X» ' VI  =  SF4.1) 

240  FORMAT <2X» 'VI  =  ' »F4 . 1 .' (MPH) '. 3X »' DRIVER  COEFFICIENTS  FOR'. 

4  '  OUTERTRACKS '»F6.1>3('»'»F6.1)»/.43X» 

4  'INNER  TRACKS'. F6. 1 >3( ' » '  »F6 . 1  > > 

245  FORMAT (28X.A76) 

246  FORMAT ( 1H1 ) 

250  F0RMAT(48X. 'TRANSIENT  TURNING  MOTION') 

C 

END 
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*  * 

*  THIS  SUBROUTINE  CALCULATES  THE  SLIP  RADII  AND  OFFSET.  * 

*  (TECHNICAL  REPORT  GL-79-6)  * 

*  * 

mtmutuumiuMtmtttmtmmmtmmtuuumm 

SUBROUTINE  AGIL(*> 

PARAMETER  (  N  =  10  »  N1  =  N+l  ) 

REAL  L»JACK»IZ»IZB 
LOGICAL  LOGREtLVX 
CHARACTER  TESTN*40»DATE*8 
SAVE 

DATA  PTA»PTAD»PTAPM» IFIRST  /0.95  »  0.95  »  0.75  »  0/ 

DIMENSION  T1(N1) »T2(N1 ) »Q1 (N1 ) »Q2(N1 > »DX(N1 > 

COMMON  /INOUT  /  EfVXI *VX2iV.RIOL»LOGREfLVX 

COMMON  /IN  /  DTIM1 »DTIN2tWPRV»UDPRVt VXPRV* VYPRV»TESTN>DATE 
COMMON  /VPARAM/  UT.L.HiDiB.SL»CX.THETAA»THETAB>IZ»VCIl 
COMMON  /SPARAM/  Cl » A. SM. SN.CtSXI tCB.CLAMDA.PHI »GfSF .ETA.CH1 
COMMON  /OUT  /  FCX»FCY»MSX1fVSX2»PTE»PTS»WD»W»VXfVYf 
l  P  fZI fZ2  f IDIR  f VXDfVYDfPT 1  fPT2 

COMMON  /MI SC  /  GRAVfDRADfVMPfRDEG 

FUNCTION  FOR  R  (INPUT  E> 

FR(Z1 fZ2fE)  =  <E*<B0L-2.*Z2>+BQL+2.*Z1>  /  (E-l.)*.5 
FUNCTION  FOR  E  (INPUT  RO) 

FE(Z1»Z2»R0L)  =  (Z1+BOL/2.+ROL)  /  (Z2-B0L/2.+R0L) 


FUNCTIONS  FOR  NORMAL  STRESSES  WITHOUT  TRACK  TENSION  TERM 

FR1B(X)  =  CZ2  +  CXCZT6*X  -  HOBFCY  -  H0BFCX*X 
FR2B(X)  =  CZ2  +  CXCZT64X  +  HOBFCY  -  H0BFCX4X 

WHERE 

HOBFCY  =  HQL/BOLK  -FCY-SZSC) 

HOBFCX  =  6*H0L*(FCX+SZCC> 

CXCZT6  =  6*CX0UC0S(ETA> 


FUNCTIONS  FOR  TRACK  TENSION  TERMS 

FSO(X)  *  AL3PRX  +  ( ALSAN1BX+AL1AN1 )*DU1  +  AN14DU4 
FSI(X)  *  AL3PRX  +  (ALSA2P*X+AL1AN2)*DU1  +  AN2A2P4DU2 
l  -  ( ALSAN24X-AL1  A2F')tDU3 
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C 

C 


EQUATIONS  FOR  TRACK  TENSION 
FA(TIS)  =  AN1  -  STETDHT1S 

GA(T2S)  =  DUM1*<AN2-STHETA*T2S>  +  DUM2*(AN2P+STETD*T2S) 
FA2(T1S)  =  (FA<T1S)-AA)  /  DT 
GA2( T2S)  =  <GA<T2S)-BB)  /  DT 


FUNCTION  FOR  TOTAL  SHEAR  STRENGTH 
EXPS(X)  =  EXP(AMAXl<X»-88. ) ) 

FSM(XCDfSIG)  =  (ABAR+SHBAR*EXPS<-SN*SIG>)*(1.-SXI)  + 
t  (CBAR+SIG*TPHI)*SXI  +  XCD 


EQUATIONS  OF  MOTION 

FB<TS)  =  TS  -  SF  -  FCX  -  SZCC 
GB(QS)  =  QS  +  FCY  t  SZSC 

HB(QXrTD)  =  QX  +  B0L/2.4TD  +  SF»HOL*<FCY+SZSC>  +  IZBtUDD 
F2<TS)  *  (FB(TS)-AV)  /  DT 
G2(QS)  =  (GB(QS)-BV)  /  DT 
H2(QXf TD)  =  (HB(QX»TD)-CV>  /  DT 


IF  ( IFIRST  .NE.  0)  GO  TO  20 
IFIRST  =  1 

DT  =  .0005 
CHI  =  CHIUDRAD 
DX ( 1 )  =  0. 

CONT  =  l.E-5 
HPF  =  kJT/6600. 

HPFE  =  0.5/PTA 
EPREV  *  E 
P  =  0. 

Z1  =  .1 
Z2  =  -.1 


VARIABLES  RELATED  TO  VEHICLE  PARAMETERS 


ALSQOU  =  L*L/«T 
HOB  =  H/B 
DOL  =  D/L 
IZB  =  IZ/L/WT 
HOL  =  H/L 
HOLT 6  =  6.IH0L 
BOL  =  B/L 
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CXOL  =  CX/L 

USED  FOR  TRACK  TENSION  ONLY  (SL  .NE.  O.) 

IF<SL)  69i73»69 

69  TETD1  =  THETADIDRAD 
THETA1  =  THETAA*DRAD 
STETD  =  SIN(TETDl) 

STHETA  =  SIN(THETAl) 

AN2  =  STHETA 

AN1  =  STETD 
AN2P  =  -STETD 
SLOL  =  SL/L 
AL2  =  SL0L*2. 

N  =  SL0L*N+.5 
Ml  =  M+l 
M2  =  Nl-M 
DX<2)  =  SLOL/M 
DO  70  I  =  3iMl 

70  DX< I )  =  DX<2) 

DX(M+2)  =  (1.-AL2)  /  <N-2*M> 

DO  71  I  =  M+3,M2 

71  DX(I)  =  DX(M+2) 

DO  72  I  =  M2+1 >N1 

72  DX< I )  =  DX<2) 

Ml  =  Ml  U 
M2  =  M2  +  1 
AL3  =  3.  -  AL2 
ALS  =  2./SL0LM2 

AL1  =  ((l.-SL0L)/SL0L)**2 
GO  TO  75 

USED  ONLY  FOR  ZERO  TRACK  TENSION  (SL  .Ed.  0.) 

73  DXT  =  l./N 

DO  74  I  =  2>N1 

74  DX< I )  =  DXT 


VARIABLES  RELATED  TO  SOIL  MODEL  PARAMETERS 


75  CONTINUE 

PHI1  =  PHI*DRAD 
TPHI  =  TAN(PHIl) 

CBAR  =  DOL*ALSQOU*C 
ABAR  =  DOL*ALSQOW*A 
SMBAR  =  BOLf ALSQ0W4SM 
GBAR  =  DOL*G*ALSQOU*L 
CDBAR  =  BOL*CD*ALSOOU 
CLMDAB  =  CLAMDAtL 
DI  =  SF 

CZ  =  COS(ETA*DRAD) 

CZ2  *  CZ/2. 

SZ  =  SIN(ETAIDRAD) 
CXCZT6  =  CX0L*CZ*6. 


AGIL0113 
AGIL0114 
AGIL01 15 
AGIL0116 
AGIL0117 
AGIL0U8 
AGIL0119 
AGIL0120 
AGIL0121 
AGIL0122 
AGIL0123 
AGIL0124 
AGIL0125 
AGIL0126 
AGIL0127 
AGIL0128 
AGIL0129 
AGIL0130 
AGIL0131 
AGIL0132 
AGIL0133 
AGIL0134 
AGIL0135 
AGIL0136 
AGIL0137 
AGIL0138 
AGIL0139 
AGIL0140 
AGIL0141 
AGIL0142 
AGIL0143 
AGIL0144 
AGIL0145 
AGIL0146 
AGIL0147 
AGIL0148 
AGIL0149 
AGIL0150 
AGIL0151 
AGIL0152 
AGIL0153 
AGIL0154 
AGIL0155 
AGIL0156 
AGIL0157 
AGIL0158 
AGIL0159 
AGIL0160 
AGIL0161 
AGIL0162 
AGIL0163 
AGIL0164 
AGIL0165 
AGIL0166 
AGIL0167 
A6IL0168 


nnn  o  r>  o  n  o  o  o  r>  o  o  o  r>  o  o  o  ooooo 


AGIL0169 

. . . . - . — .  AGIL0170 

DETERMINE  DIRECTION  VEHICLE  IS  TURNING  AGIL0171 

AGIL0172 
AGIL0173 
AGIL0174 
AGIL0175 
AGIL0176 
AGIL0177 
AGIL0178 
AGIL0179 
AGIL0180 
AGIL0181 
AGIL0182 
AGIL0183 
AGIL0184 
AGIL0185 
AGIL0186 
AGIL0187 
AGIL0188 
AGIL0189 
AGIL0190 
AGIL0191 
AGIL0192 
AGIL0193 
AGIL0194 
AGIL0195 
AGIL0196 
AGIL0197 
AGIL0198 
AGIL0199 
AGIL0200 
AGIL0201 


ITERATE  TO  FIND  SOLUTIONS  AGIL0202 

— - AGIL0203 

AGIL0204 

DO  100  IL  =  1 >36  AGIL0205 

AGIL0206 

K  =  1  AGIL0207 

8  PSO  =  P**2  AGIL0208 

IF  (LOGRE)  GO  TO  510  AGIL0209 

AGIL0210 

FOR  CASES  WHERE  E  IS  INPUT  AGIL02U 

AGIL0212 

ROL  =  FR(Z1»Z2?£)  AGIL0213 

RIOL  =  SIGN(SGRT(PSQ+R0L**2)»R0L)  AGIL0214 

GO  TO  520  AGIL0215 

AGIL0216 

FOR  CASES  WHERE  RI  IS  INPUT  <RI=RO  FOR  STEADY-STATE)  AGIL0217 

AGIL0218 

510  CONTINUE  AGIL0219 

RIOLSQ  =  RI0L4RI0L  AGIL0220 

IF  < F'SQ  .GE.  RIOLSQ)  RETURN  1  AGIL0221 

ROL  =  SIGN(SQRT(RIOLSQ-PSQ)»RIOL)  AGIL0222 

E  =  FE(Z1»Z2.R0L>  AGIL0223 

C  AGIL0224 


20  CONTINUE 
IDIR  =  1 

IF  (LOGRE)  GO  TO  47 
IF(E-1.)  41»44*47 

41  IF(E)  42»43f43 

42  E  =  -1,/E 

IF  (E-l . )  43»44»47 

43  CONTINUE 
IDIR  =  -1 
GO  TO  47 

44  CONTINUE 

IF (EPREV-1 , )  45» 46, 46 

45  E  =  0.9999 

IF  (LVX)  VX1  =  VX2*E 
IDIR  =  -1 
GO  TO  47 

46  E  =  1.0001 

IF  (LVX)  VX2  =  VX1/E 

47  CONTINUE 
EPREV  =  E 


nnn 


520 

CONTINUE 

AGIL0225 

IF  (LOX)  GO  TO  530 

AGIL0226 

WDL  =  O/RIOL 

AGIL0227 

0X2  =  WDL*<B0L+Z1-Z2)/<E-1.) 

AGIL0228 

0X1  =  0X24E 

AGIL0229 

GO  TO  540 

AGIL0230 

AGIL0231 

530 

CONTINUE 

AGIL0232 

WDL  =  0X2KE-1  • )  /  (B0L+Z1-Z2) 

AGIL0233 

AGIL0234 

o 

*r 

m 

WDL02  =  WDL*. 5 

AGIL0235 

UD  =  WDL/L 

AGIL0236 

OX  =  R0L4UDL 

AGIL0237 

Or  =  PtUDL 

AGIL0238 

UDD  =  (WD-UDPR0XDTIM1 

AGIL0239 

OXD  =  (OX-OXPRO)*DTIhl 

AGIL0240 

OYD  =  (OY-OYPROXBTIMl 

AGIL0241 

FCY  =  <OX#UO-UYD)/GRAU 

AGIL0242 

FCX  =  <UY*UD+OXD)/GRAU 

AGIL0243 

U  =  WF'RO  +  ( WD+WDPRO )  *DT  IH2 

AGIL0244 

CC  =  COS(CHI+W) 

AGIL0245 

SC  =  SIN(CHItW) 

AGIL0246 

SZSC  =  SZ*SC 

AGIL0247 

SZCC  =  SZ*CC 

AGIL0248 

HOBFCY  =  HOBt(-FCY-SZSC) 

AGIL0249 

HOBFCX  =  H0LT6*(FCX+SZCC) 

AGIL0250 

Z1SQ  =  Z1*Z1 

AGIL0251 

Z2SQ  =  Z2*Z2 

AGIL0252 

IF (0)  77»78i77 

AGIL0253 

77 

TEH1  =  WDL02/OX1 

AGIL0254 

TEM2  =  WDL02/0X2 

AGIL0255 

GO  TO  79 

AGIL0256 

AGIL0257 

78 

TEH1  =  0, 

AGIL0258 

TEH2  =  0. 

AGIL0259 

AGIL0260 

79 

CONTINUE 

AGIL0261 

0SX2  =  Z2*WDL 

AGIL0262 

K2  =  4 

AGIL0263 

AGIL0264 

USED  FOR  TRACK  TENSION  ONLY  (SL  .NE.  0.) 

AGIL0265 

AGIL0266 

IF(SL)  80i85t 80 

AGIL0267 

80 

IF <0SX2)  81 i82t 82 

AGIL0268 

81 

AN2  =  0. 

AGIL0269 

DUhl  =  0. 

AGIL0270 

DUM2  =  1. 

AGIL0271 

GO  TO  83 

AGIL0272 

82 

AN2P  *  0. 

AGIL0273 

DUM1  *  1. 

AGIL0274 

DUN2  *  0. 

AGIL0275 

83 

CONTINUE 

AGIL0276 

K2  =  1 

AGIL0277 

IL2  =  1 

AGIL0278 

84 

ALSAN1  =  ALStANl 

AGIL0279 

AL1AN1  =  ALUAN1 

AGIL0280 

o  n  o  o  r> 


AL3*<AN2-AN2P-AN1> 
ALS*AN2P 
AL14AN2P  4  AN2 
AN24AN2P 
ALS*AN2 

AL14AN2  +  AN2P 


AL3PR  = 

ALSA2P 

AL1AN2 

AN2A2P 

ALSAN2 

AL1A2P 

DU1=1. 

DU2=0. 

DU3=0. 

DU4=0. 


85  CONTINUE 
X  =  -.5 

UP1  =  .5  -  P  -  CXOL 

UP1S  =  UP1**2 

UP1Z1S  =  SORT (UP1S4Z1SG) 

UP1Z2S  =  SQRT(UP1S4Z2SQ) 

ARG1  =  AMAX1 (UP14UP1Z1S , 1 .E-20 ) 

ARG2  =  ANAX1 (UP14UP1Z2S, 1 < E-20) 

DIBP  =  UP1*UP1Z1S  +  Z1SQ*AL0G(ARG1) 
D2BP  =  UP14UP1Z2S  4  Z2SQIAL0G ( ARG2 ) 


INTEGRATE  ALONG  TRACKS  FOR  COMPONENTS  OF  SHEAR  STRESS 


DO  10  I  =  1 ,N1 

X  =  X  4  DX  Cl) 

RIB  =  FR1B(X) 

R2B  =  FR2B(X) 

USED  FOR  TRACK  TENSION  ONLY  <SL  ,NE.  0.) 

IF(SL)  32.37.32 
IF(I-Ml)  36.33.34 
DU1=0. 

DU2=1 . 

DU4*1 . 

IF ( I-M2)  36.35,36 
DU2=0. 

DU3=1. 

CONTINUE 

AL3PRX  =  AL3PR4X 
RIB  =  RIB  4  FSO(X) 

R2B  =  R2B  4  FSI(X) 

CONTINUE 

XNP  =  X  -  P  -  CXOL 

XMPSQ  =  XMP**2 

RTZ1  =  SORT (XMPSQ4Z1S0) 

RTZ2  =  SORT ( XMPSQ4Z2SQ) 

SGI  =  XMP/RTZ1 
CGI  =  Z1/RTZ1 
SG2  =  XMP/RTZ2 
CG2  =  Z2/RTZ2 


AGIL0281 

AGIL0282 

AGIL0283 

AGIL0284 

AGIL0285 

AGIL0286 

AGIL0287 

AGIL0288 

AGIL0289 

AGIL0290 

AGIL0291 

AGIL0292 

AGIL0293 

AGIL0294 

AGIL0295 

AGIL0296 

AGIL0297 

AGIL0298 

AGIL0299 

AGIL0300 

AGIL0301 

AGIL0302 

AGIL0303 

AGIL0304 

AGIL0305 

AGIL0306 

AGIL0307 

AGIL0308 

AGIL0309 

AGIL0310 

AGIL0311 

AGIL0312 

AGIL0313 

AGIL0314 

AGIL0315 

AGIL0316 

AGIL0317 

AGIL0318 

AGIL0319 

AGIL0320 

AGIL0321 

AGIL0322 

AGIL0323 

AGIL0324 

AGIL0325 

AGIL0326 

AGIL0327 

AGIL0328 

AGIL0329 

AGIL0330 

AGIL0331 

AGIL0332 

AGIL0333 

AGIL0334 

AGIL0335 

AGIL0336 


fcsr 
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DD1B  =  WD*RTZ1 
DD2B  =  WD*RTZ2 

ARG1  =  AMAX1  ( XMP+RTZ1  »1.E— 20) 

ARG2  =  AMAX1(XHP+RTZ2>1 .  E-20) 

DIB  =  DI  +  TEM1*<D1BP-XMP*RTZ1-Z1SQ*AL0G(ARG1>) 
D2B  =  DI  +  TEM2*(D2BP-XMP*RTZ2-Z2SQ*AL0G(ARG2)> 
GD1B  =  DIBtGBAR 
GD2B  =  D2B*GBAR 

CCD1  =  CDBAR  *  <1.0-EXPS(-CLMDAB*ABS(DB1B>)> 
CCD2  =  CDBAR  *  ( 1 ,0-EXPS<-CLHDAB*ABS(DD2B> ) ) 
TAU1  =  FSM(CCD1»AMAX1(R1B»0. )) 

TAU2  =  FSM(CCD2»AHAX1(R2B»0. ) ) 

TQ1  =  GD1B*TAU1  /  <TAU1+ABS(GD1B> > 

TQ2  =  GD2B*TAU2  /  ( TAU2+ABS ( GD2B ) ) 

T 1  <  I )  =  TQ 1 *CG  J 
T2 ( I )  =  TG2*CG2 
Q 1  ( I )  =  TQUSG1  +  TQ2ISG2 
02(1)  =  01 (I)*(X-CXOL) 

10  CONTINUE 


T1S  =  0. 

T2S  =  0. 

DO  21  I  =  2»N1 
DX02  =  DX(I)*.5 

T1S  =  T IS  +  (T1(I)+T1(I-1))*DX02 
21  T2S  =  T2S  +  <T2(I)+T2(I-1))*DX02 
T1S  =  T1S*IBIR 
T2S  =  T2S*IDIR 
GO  TO  (87*88.89.90)  ,K2 

USED  FOR  TRACK  TENSION  ONLY  (SL  .NE.  0.) 

87  AA  =  FA(TIS) 

BB  =  GA(T2S) 

IF(ABS(AA) .LT.CONT  .AND.  ABS(BB)  .LT.CONT)  GO  TO  90 
K2  =  2 

AN2  =  AN2  +  DUN1IDT 
AN2P  =  AN2P  +  DUH2*DT 
GO  TO  84 

88  FAX  =  FA2(T1S) 

GAX  =  GA2(T2S) 

AN2  =  AN2  -  DUMUDT 
AN2P  =  AN2P  -  DUM2CDT 
AN1  =  AN1  +  DT 
K2  =  3 
GO  TO  84 

89  FAY  =  FA2(T1S) 

GAY  =  GA2(T2S) 

AN1  =  AN1  -  DT 

RJAC2  =  FAX*GAY  -  FAY4GAX 
IF< ABS(RJAC2)  .LT.  l.E-37)  RETURN  1 
DELTTT  =  ( -AA*GAY+BB*FAY )  /  RJAC2 
DELTT3  =  (-BB*FAX+AA*GAX)  /  RJAC2 


AGIL0337 

AGIL0338 

AGIL0339 

AGIL0340 

AGIL0341 

AGIL0342 

AGIL0343 

AGIL0344 

AGIL0345 

AGIL0346 

AGIL0347 

AGIL0348 

AGIL0349 

AGIL0350 

AGIL0351 

AGIL0352 

AGIL0353 

AGIL0354 

AGIL0355 

AGIL0356 

AGIL0357 

AGIL0358 

AGIL0359 

AGIL0360 

AGIL0361 

AGIL0362 

AGIL0363 

AGIL0364 

AGIL0365 

AGIL0366 

AGIL0367 

AGIL0368 

AGIL0369 

AGIL0370 

AGIL0371 

AGIL0372 

AGIL0373 

AGIL0374 

AGIL0375 

AGIL0376 

AGIL0377 

AGIL0378 

AGIL0379 

AGIL0380 

AGIL0381 

AGIL0382 

AGIL0383 

AGIL0384 

AGIL038S 

AGIL0386 

AGIL0387 

AGIL0388 

AGIL0389 

AGIL0390 

AGIL0391 

AGIL0392 


AN2  =  AN2  +  BUMUBELTTT 
AN2P  =  AN2P  +  DUM2*DELTTT 
AN1  =  AN1  +  DELTT3 
IL2  =  IL2  +  1 
K2  =  1 

IF  (IL2  .LT.  36)  GO  TO  84 
RETURN  1 

90  CONTINUE 
QS  =  0. 

QX  =  0. 

DO  23  I  =  2.N1 
DX02  =  DX< I  )*. 5 

QS  =  QS  +  (Q1<I)+Q1(I-1))*DX02 
23  QX  =  QX  +  (Q2( I >+Q2< 1-1 ) ) *DX02 
QS  =  QS*IDIR 
QX  =  QX  *  I  III  R 
TS  =  T1S  +  T2S 
TD  =  T2S  -  T1S 
GO  TO  (30,40.50,60)  ,K 

30  AV  =  FB(TS) 

BV  =  GB(QS) 

CV  =  HB(QX.TD) 

IF(ABS(AV) .LT.CONT  .AND.  ABS(BV) 
X  .AND.  ABS(CV) 

P  =  P+DT 
K  =  2 
GO  TO  8 

40  FX  =  F2(TS> 

GX  =  G2<  QS ) 

HX  =  H2(QX»TD) 

P  «  P  -  DT 
Z1  =  Z1  +  DT 
K  =  3 
GO  TO  8 

50  FY  =  F2(TS) 

GY  =  G2(QS) 

HY  =  H2(QX,TD) 

Z1  =  Z1  -  DT 
Z2  “  Z2  +  DT 
K  =  4 
GO  TO  3 

60  FZ  =  F 2 ( T S ) 

GZ  =  G2<  QS  > 

HZ  =  H2(QX,TD) 

Z2  =  Z2  -  DT 
A1F  =  GY*HZ  -  GZ*HY 

B1F  =  FY*HZ  -  FZtHY 

C1F  =  FY*GZ  -  FZ*GY 

A2F  =  BV*HX  -  CV*GX 

C2F  =  AV*HX  -  CV*FX 

C2F  =  AV*GX  -  BV*FX 


AGIL0393 
AGIL0394 
AGIL0395 
AGIL0396 
AGIL0397 
AGIL0398 
AGIL0399 
AGIL0400 
AGIL040I 
AGIL0402 
AGIL0403 
AGIL0404 
AGIL0405 
AGIL0406 
AGIL0407 
AGIL0408 
AGIL0409 
AGIL0410 
AGIL041 1 
AGIL0412 
AGIL0413 
AGIL0414 
AGIL0415 
AGIL0416 

LT.CONT  AGIL0417 

LT.CONT)  GO  TO  120  AGIL0418 

AGIL0419 

AGIL0420 

AGIL0421 

AGIL0422 

AGIL0423 

AGIL0424 

AGIL0425 

AGIL0426 

AGIL0427 

AGIL0428 

AGIL0429 

AGIL0430 

AGIL0431 

AGIL0432 

AGIL0433 

AGIL0434 

AGIL0435 

AGIL0436 

AGIL0437 

AGIL0438 

AGIL0439 

AGIL0440 

AGIL0441 

AGIL0442 

AGIL0443 

AGIL0444 

AGIL0445 

AGIL0446 

AGIL0447 

AGIL0448 
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OOOOOOOOO 


c 

JACK  *  FXtAlF  -  GXBB1F  4  HXBC1F 
IF(ABS( JACK) <LT < 1 <E*37)  RETURN  1 
C 

P  *  P  ♦  (BVBB1F  -  AVBA1F  -  CV*C1F)/JACK 
21  8  ZH  (FZBA2F  -  GZ*B2F  +  HZ*C2F)/JACK 
Z2  =  Z2  +  <GY*B2F  -  FY*A2F  -  HY*C2F)/JACK 
C 

100  CONTINUE 
RETURN  1 

C  — . - . . . 

C 

120  CONTINUE 

PT1  *  VXIBTlStHPF 
PT2  =  VX2*T2S*HPF 
PTS  =  PT1  +  PT2 
IF  <  E-l . )  140.130.130 
C 

130  PTE  =  HPFE*(PT1/E*<E+1.+<E-1.)/PTAPMH 
I  PT2*(E+1 .  -PTADXE-1  .  )/PTAPM) ) 

GO  TO  150 
C 

140  ONEOE  *  l./E 

PTE  =  HPFE*<PT2*E#<0NE0E41.+<0NE0E-1,)/PTAPH>4 
X  PT 1 * ( ONEOE  +  1 . -PTAD* ( ONEOE-1 , ) /PTAPH ) ) 

C 

150  CONTINUE 

VSX1  =  Z1XUDL 
C 

RETURN 
END 

************************************************************* 
*  * 

*  THIS  SUBROUTINE  CALCULATES  THE  VEHICLE  CONE  INDEX  * 

*  FOR  ONE-PASS  TRAFFIC  USING  THE  MOBILITY  INDEX.  * 

*  < TECHNICAL  REPORT  GL.-79-6)  * 

*  * 
umttummmmtmmmmmmmmmtmmmm 

SUBROUTINE  MIVCI 
C 

REAL  L 

COMMON  /VPARAM/  UT.L.H.D.B.SL.CX.THETAA.THETAD. I Z » VC 1 1 
COMMON  /MISC  /  GRAV.DRAD. VMP.RDEG 
C 

PRINT  *» 

PRINT  *»  'FOR  TRNT »  0  =  MANUAL.  1  =  AUTOMATIC' 

PRINT  *.  'INPUT  GH.NB.TSL.HP. TRNT' 

READ  *.  GH.NB.TSL.HP. TRNT 
PRINT  *. 

C 

CPF  =  HT/(D*L) 

C 

WF  =  1.0 

IFIWT.GE. 50000.0  .AND.  WT.LT. 70000.0)  WF  =  1.2 


AGIL0449 
AGIL0450 
AGIL0451 
AGIL0452 
AGIL0453 
A6IL0454 
AGIL0455 
AGIL0456 
AGIL0457 
AGIL045B 
AGIL0459 
AGIL0460 
AGIL0461 
AGIL0462 
AGIL0463 
AGIL0464 
AGIL0465 
AGIL0466 
AGIL0467 
AGIL0468 
AGIL0469 
AGIL0470 
AGIL0471 
AGIL0472 
AGIL0473 
AGIL0474 
AGIL0475 
AGIL0476 
AGIL0477 
AGIL0478 
AGIL0479 
MIVCI001 
MIVCI002 
MIVCI003 
MIVCI004 
MIVCI005 
MI VCI006 
MIVCI007 
MIVCI008 
MIVCI009 
MIVCI010 
MIVCI01 1 
MIVCI012 
MIVCI013 
MIVCI014 
MIVCI015 
MIVCI016 
MIVCI01 7 
MIVCI018 
MIVCI019 
MIVCI020 
MIVCI021 
MIVCI022 
MIVCI023 
MIVCI024 
MIVCI025 
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IFIWT.GE. 70000,0  .AND.  WT.LT. 100000.0)  WF  =  1.4 
IF<WT  .GE.  100000.0)  WF  =  1.8 

TRKF  =  D/100. 

GF  =  1.0 

IF <GH  .GT.  1.5)  GF  =  1.1 

BF  =  <WT/10.)/<NB*TSL*D) 

CF  =  1.5*L/<NB-1)/TAN(THETAA*DRAD)/10. 

HPT  =  HP/(WT/2000. ) 

EF  =  1.0 

IF(HPT  .LT.  10.)  EF  =  1.05 
TRNF  =1.05 

IF (TRNT  .GT.  0.)  TRNF  =  1.0 

AMI  =  (<CPF*UF)/(TRKF*GF)+BF-CF)  *  EF  *  TRNF 


VCI1  =  7.0  ♦  0.24AMI  -  39.2/<AMI+5. 6) 


RETURN 

END 


* 

THIS  SUBROUTINE  CALCULATES  THE  MOBILITY  CONE  INDEX.  * 

TO  SIMULATE  THE  STANDARD  WES  CONE  INPUT  * 

CL  =  1.48  (IN)  »  DI  =  0.799  (IN)  * 

(MISCELLANEOUS  PAPER  SL-81-4)  * 

* 


M*tt*tt***ttttt***tmt*m*tmt*m*mm*************tm 

SUBROUTINE  CONE 

DOUBLE  PRECISION  TL»AG. AGM.AA,A1»A2>A3>A4»A5,A6>A12 
COMMON  /SPARAM/  CI»A»SM»SN»C»SXIiCD»CLAMDA»PHI»G»SF»ETA»CHI 
COMMON  /MISC  /  GRAV » DRAB  * VMP » R  BEG 

PRINT  *> 

PRINT  *»  'TO  SIMULATE  STANDARD  WES  CONEf  INPUT  CL= 1 . 48 r DI=0 . 79 ' 
PRINT  *»  'FOR  TYPICAL  GAMA.Z,  INPUT  GAMA=0.07»Z=6. ' 

PRINT  *t  'INPUT  CLfDI.GAMAjZ' 

READ  *,  CLfDIfGAMAfZ 

TL  =  CL 
AG  =  4*G 
AGM  =  GAMA 
TAl  =  DI/2./CL 
PHI 1  =  PHItDRAD 
SPHI  =  SIN(PHIl) 


MIVCI026 
MIVCI027 
MIVCI028 
MIVCI029 
MIVCI030 
MIVCI031 
MIVCI032 
MIVCI033 
MIVCI034 
MI  VCI035 
MIVCI036 
MIVCI037 
MIVCI036 
MIVCI039 
MIVCI040 
MIVCI041 
MIVCI042 
MIVCI043 
MI VCI044 
MIVCI045 
MIVCI046 
MIVCI047 
MIVCI048 
MIVCI049 
MIVCI050 
MIVCI051 
MIVCI052 
MIVCI053 
C0NE0001 
C0NE0002 
C0NE0003 
CONE0004 
C0NE0005 
C0NE0006 
C0NE0007 
CONE0008 
C0NE0009 
C0NE0010 
C0NE0011 
C0NE0012 
C0NE0013 
C0NE0014 
C0NE0015 
CONE0016 
C0NE0017 
C0NE0018 
CONE0019 
C0NE0020 
CONE0021 
C0NE0022 
C0NE0023 
C0NE0024 
C0NE0025 
C0NE0026 
C0NE0027 
CONE0028 


o  o  o  non  onooo  oor> 


C 

C 


TPHI  =  TAN<PHI1 ) 

ELP  =  4.<SPHI/<1.4SPHI)/3. 
ELP2  =  2.  -  ELP 
ELP3  =3.  -  ELP 

IF(PHI)  30»20»30 

20  CONTINUE 


Cl  =  4.<C/3.<<1.+DL0G<AG/C>)  + 

t  2.<C<TL/DI  +  AGN<<<Z+TL)-2.<TL/3.) 


RETURN 


30  CONTINUE 

A1  =  AGM*TPHI 
AA  =  A1<<Z4TL)  +  C 
A12  =  Al«2 

A2  *  A12<TPHI<< DI/2. X<2 
A3  =  2«<TAL<< 1 . 4SPHI ) <AG<<ELP 
A4  =  3.<<TAL4TPHI)/(3.-SPHI) 

A5  =  AA«ELP3 

A6  =  <AA4ELP2<A1<TLXDNAX1  (0 . » ( AA-TL<AGN<TPHI )  )«ELP2 


Cl  =  -C/TPHI  +  A37A2<A4<<  AS- A6 ) /ELP2/ELP3 


RETURN 

END 


CQNE0029 

CQNE0030 

C0NE0031 

C0NE0032 

C0NE0033 

C0NE0034 

C0NE0035 

C0NE0036 

C0NE0037 

C0NE0038 

C0ME0039 

C0NE004O 

C0NE0041 

C0NE0042 

C0NE0O43 

C0NE0044 

C0NE0045 

C0NE0046 

CONE0047 

C0NE0048 

C0NE0049 

C0NE0050 

CONE0051 

CONE0052 

C0NE0053 

C0NE0054 

C0NE0055 

CONE0056 

C0NE0O57 

C0NE00S8 

C0NE0059 

CQNE0060 

C0NE0061 

CONE0062 

C0NE0063 

C0NE0064 

C0NE0065 
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APPENDIX  B:  GLOSSARY 


B.1  MAIN  PROGRAM 

A 

B 

C 

CD 

CHI 

Cl 

CLAMDA 

CONE 

CX 

D 

DATE 

DATIM 

DRAD 

ETA 

FPRAM 


FXOPT 

G 

GRAV 

H 

IRUN 

IZ 


L 

MIVCI 


Material  constant  in  failure  envelope,  psi. 

Track  tread  (distance  between  center  lines  of  tracks), 
in. 

Static  cohesive  component  of  shear  strength,  psi. 

Added  cohesive  strength  due  to  dynamic  loading,  psi. 
Directional  angle  of  the  vehicle  on  sloping  terrain,  deg. 
WES  cone  index. 

Material  constant  related  to  rate  effect,  sec/in. 
Subroutine  to  calculate  Cl. 

Abscissa  of  the  center  of  gravity  of  the  vehicle,  in. 
Track  width,  in. 

Character  variable  for  the  current  date. 

System  subroutine  for  obtaining  the  current  date  and 
time. 

Factor  for  converting  from  degrees  to  radians. 

Angle  of  sloping  terrain,  deg. 

System  subroutine  for  setting  the  reflexive  read 
characters  that  are  sent  to  a  terminal  for  requesting 
input  and  the  line  length  for  formatted  output  directed 
to  a  terminal. 

System  subroutine  for  monitoring  FORTRAN  execution 
errors. 

Shear  modulus,  psi/in. 

2 

Acceleration  due  to  gravity,  in/sec  . 

Height  of  the  center  of  gravity,  in. 

Variable  used  to  determine  the  type  of  run  (steady-state 
or  transient  motion). 

Mass  moment  of  inertia  of  the  vehicle  about  an  axis 
passing  through  its  center  of  gravity  and  parallel  to  the 

p 

Z  axis,  in-lb/sec  . 

Length  of  track  in  contact  with  the  ground,  in. 

Subroutine  to  calculate  VCI1  (vehicle  cone  index  for  one 
pass) . 
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PTIME 


=  Octal  constant  used  by  subroutine  FPARAM. 

=  Parameter  variable  for  setting  dimension  of  arrays. 

=  Angle  of  internal  friction,  deg. 

=  Processor  time  at  start  of  calculation,  hr. 

=»  System  subroutine  for  obtaining  the  central  processor 


PTU  =  Processor  time  used,  min. 

RDEG  =  Factor  for  converting  from  radians  to  degrees. 

SF  =  Coefficient  of  rolling  resistance. 

SL  =  Distance  between  center  lines  of  adjacent  wheels,  in. 
SM  =  Material  constant  in  failure  envelope,  psi. 

SN  =  Material  constant  in  failure  envelope,  1/psi. 

STEADY  =  Subroutine  for  steady-state  turning  motion. 


TESIN  = 
THETAA  = 
THETAD  = 


Variable  used  to  determine  the  type  of  soil  model  used. 
Test  title  of  up  to  40  characters. 

Approach  angle  of  the  track  envelope,  deg. 

Departure  angle  of  the  track  envelope,  deg. 


TIM  =  Current  time,  used  only  for  subroutine  DATIM. 

TRANS  =  Subroutine  for  transient  motion. 

VCI1  -  Vehicle  cone  index  for  one  pass. 

VMP  =  Factor  for  converting  from  inches  per  second  to  miles  per 
hour. 

WT  =  Weight  of  the  vehicle,  lb. 


B.2  SUBROUTINE  STEADY 


DTIM1 


Subroutine  for  solving  the  equations  of  motion. 

Steering  ratio  increment. 

Turning  radius  decrement,  ft. 

Set  to  zero  for  steady-state  turning  motion. 

Velocity  increment,  mph  (converted  internally  to  in/sec.) 
Steering  ratio. 

E  at  the  power  limit. 

Array  for  storing  E. 

E  at  the  stability  limit. 

Forward  acceleration,  g's. 

FCX  at  the  power  limit,  g's. 


kS 
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Array  for  storing  FCX. 


=  FCX  at  the  stability  limit,  g's. 


Lateral  acceleration,  g's. 
FCY  at  the  power  limit,  g's. 
Array  for  storing  FCY. 


=  FCY  at  the  stability  limit,  g's. 


ICASE  = 


Horsepower  available,  hp. 
Case  number:  1,  2,  3,  4  or  ? 


=  Flag  used  for  power  limit  check. 


=  Loop  counter. 


INPUT  =  Character  variable  for  keyboard  input. 


LOGRE  = 


Length  of  track  in  contact  with  the  ground,  in. 
Character  variable  for  type  of  stability  limit. 
Character  array  for  the  case  menu. 

Switch  variable;  set  false  if  E  is  input,  true  if  R0  is 


input. 


=  Switch  variable  for  stability  check;  set  true  when 


stability  limit  is  reached. 


Set  false  for  steady-state  turning  motion. 

Parameter  variable  for  setting  dimension  of  arrays. 
Point  number  at  the  stability  limit. 


»  Maximum  number  of  points  which  can  be  calculated  for 


Cases  1  and  2. 


=  Number  of  points  to  be  calculated. 


=  Offset. 


Factor  for  interpolation. 


P  at  the  power  limit. 


Factor  for  interpolation. 


=  Array  for  storing  P. 


=  P  at  the  stability  limit. 


Power  required  by  the  sprocket  of  the  outer  track,  hp. 


=  Array  for  storing  PT1 . 


Power  required  by  the  sprocket  of  the  inner  track,  hp. 
Array  for  storing  PT2. 

Power  required  from  the  engine,  hp. 

PTE  at  the  power  limit,  hp. 
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PTEP  =  Array  for  storing  PTE. 

PTES  -  PTE  at  the  stability  limit,  hp. 

PTS  -  Total  power  required  by  the  sprockets,  hp. 

PTSF  »  PTS  at  the  power  limit,  hp. 

PTSP  =  Array  for  storing  PTS. 

PTSS  -  PTS  at  the  stability  limit,  hp. 

R0  »  Radius  of  the  trajectory  (turning  radius)  of  the  center 
of  gravity  of  the  vehicle,  ft  (converted  internally  to 
in). 

R0P  =  Array  for  storing  R0. 

R0PF  =  R0  at  the  power  limit,  ft. 

R0PS  =  R0  at  the  stability  limit,  ft. 

RDEG  *  Factor  for  converting  from  radians  to  degrees. 

RI  =  Instantaneous  radius  of  curvature,  in. 

RIOL  =  RI/L. 

V  =  Velocity  of  the  vehicle,  mph  (converted  internally  to 
in/sec) . 

VFMT  =  Character  variable  containing  format  for  input  variables. 
VMP  =  Factor  for  converting  from  in/3ec  to  mph. 

VP  =  Array  for  storing  V. 

VPF  =  Vp  at  the  power  limit,  mph. 

VPS  =  Vp  at  the  stability  limit,  mph. 

VSX1  =  Longitudinal  component  of  slip  velocity  of  the  outer 
track,  in/sec. 

VSX1 F  =  VSX1  at  the  power  limit,  mph. 

VSX1P  =  Array  for  storing  VSX1 . 

VSX1S  =  VSX1  at  the  stability  limit,  mph. 

VSX2  =  Longitudinal  component  of  slip  velocity  of  the  inner 
track,  in/sec. 

VSX2F  =  VSX2  at  the  power  limit,  mph. 

VSX2P  =  Array  for  storing  VSX2. 

VSX2S  -  VSX2  at  the  stability  limit,  mph. 

VX1  =  Longitudinal  component  of  velocity  of  the  outer  track, 
in/sec. 

VX1F  -  VX1  at  the  power  limit,  mph. 

VX1P  -  Array  for  storing  VX1. 


VX1S  =  VX1  at  the  stability  limit,  mph. 

VX2  =  Longitudinal  component  of  velocity  of  the  inner  track, 
in/sec. 

VX2F  =  VX2  at  the  power  limit,  mph. 

VX2P  =  Array  for  storing  VX2. 

VX2S  =  VX2  at  the  stability  limit,  mph. 

VXPRV,  VYPRV  =  Set  to  zero  for  steady-state  turning  motion. 

WD  =  Yaw  rate,  rad/sec. 

WDF  =  WD  at  the  power  limit,  deg/sec. 

WDP  =  Array  for  storing  WD. 

WPPRV  =  Set  to  zero  for  steady-state  turning  motion. 

WDS  =  WD  at  the  stability  limit,  deg/sec. 

Z1  =  Slip  radius  of  the  outer  track. 

ZIP  =  Previous  value  of  Z1 . 

Z2  =  Slip  radius  of  the  inner  track. 

Z2P  *  Previous  value  of  Z2. 

B.3  SUBROUTINE  TRANS 

AGIL  -  Subroutine  for  solving  the  equations  of  motion. 

ALP HD  =  Rate  of  change  in  the  side-slip  angle,  rad/sec. 

DA  »  Velocity  increment  for  power  control,  in/sec. 

DT  =  Time  increment,  sec. 

DTIM1  =  1 /DT. 

DTIM2  *=  DT/2. 

DT012  =  DTIM1/12. 

DV  =  Velocity  increment  for  stability  control,  in/sec. 

DXT  =  V  *  cos  THETA. 

DXTP  =  Array  for  storing  DXTPS. 

DXTPRV  =  Previous  value  of  DXT. 

DXTPS  »  Abscissa  of  trajectory,  ft. 

DYT  =  V  *  sin  THETA. 

DYTP  =  Array  for  storing  DYTPS. 

DYTPRV  =  Previous  value  of  DYT. 

DYTPS  =  Ordinate  of  trajectory,  ft. 

E  =  Steering  ratio. 

El  ,  E2,  E3,  E1*  =  Coefficients  for  calculating  the  steering  ratio  E. 
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El  =  Current  value  of  E. 

EP  =  Array  for  storing  E. 

FCX  =  Forward  acceleration,  g's. 

FCXP  =  Array  for  storing  FCX. 

FCY  =  Lateral  acceleration,  g's. 

FCYP  =  Array  for  storing  FCY. 

FT  =  Driver  function  used  for  sample  runs. 

HPAV  =  Horsepower  available,  hp. 

ICASE  »  Case  number. 

IDIR  =  Turn  direction  indicator  (1  for  left  turn,  -1  for  right 
turn) . 

II  =  Loop  counter. 

INPUT  =  Character  variable  for  keyboard  input. 

IPRINT  *  Print  skip  increment. 

L  «  Length  of  track  in  contact  with  the  ground,  in. 

LFMT  -  Character  array  for  the  case  menu. 

LOGRE  =  Switch  variable  (set  true  if  RI  is  input). 

LOGV  »  Switch  variable  for  velocity  check. 

LSTB  »  Switch  variable  for  stability  check. 

LVX  -  Switch  variable  (set  true  if  track  velocities  are  input) 
NP  -  Number  of  points  t>  r.  calculated. 

N2  »  Parameter  variable  for  setting  dimension  of  arrays. 

P  -  Offset. 

PP  =  Array  for  storing  P. 

PT1  =  Power  required  by  the  sprocket  of  the  outer  track,  hp. 
PT1P  =  Array  for  storing  PT1 . 

PT2  =  Power  required  by  the  sprocket  of  the  inner  track,  hp 
PT2P  =  Arr  >y  for  storing  PT2. 

PTE  «  Power  required  from  the  engine,  hp. 

PTEP  *  Array  for  storing  PTE. 

PTS  *  Total  power  required  by  the  sprockets,  hp. 

PTSP  »  Array  for  storing  PTS. 

RATIO  »  Ratio  of  power  required  to  horsepower  available. 

R0P  «  Array  for  storing  the  radius  of  trajectory  (turning 
radius)  of  the  center  of  gravity  of  the  vehicle. 

»  Factor  for  converting  from  radians  to  degrees. 


RDEG 


RI1 ,  RI2,  RI3,  =  Coefficients  for  calculating  the  instantaneous  radius  of 
RI4  curvature  (RI). 

RIOL  =  RI/L 

T  =  Time,  sec. 

THED  =  Rate  of  change  in  the  directional  angle,  rad/ sec. 

THETA  =  Directional  angle,  rad. 

THDPRV  =  Previous  value  of  THED. 

TP  =  Array  for  storing  T. 

V  =  Velocity  of  the  vehicle,  mph  (converted  internally  to 
in/ sec ) . 

VI,  V2,  V3,  V4  =  Coefficients  for  calculating  the  velocity  of  the  vehicle 

(V). 

VFMT  =  Character  variable  containing  format  for  input  variables. 
VI  =  Initial  velocity  of  the  vehicle,  mph. 

VMP  =  Factor  for  converting  from  in/ sec  to  mph. 

VP  =  Array  for  storing  V. 

VPREV  =  Previous  value  of  V. 

VSX1  =  Longitudinal  component  of  slip  velocity  of  the  outer 
track,  in/ sec. 

VSX1P  =  Array  for  storing  VSX1 . 

VSX2  =  Longitudinal  component  of  slip  velocity  of  the  inner 
track,  in/ sec. 

VSX2P  =  Array  for  storing  VSX2. 

VX  =  Longitudinal  component  of  velocity  of  the  vehicle, 
in/ sec . 

VX1  =  Longitudinal  component  of  velocity  of  the  outer  track, 
in/ sec . 

VX1 1 ,  VX12,  =  Coefficients  for  computing  the  longitudinal  component  of 

VX13,  VX1 4  velocity  of  the  outer  track  (VX1). 

VX1P  =  Array  for  storing  VX1. 

VX 2  =  Longitudinal  component  of  velocity  of  the  inner  track, 

in/sec. 

VX21 ,  VX22,  =  Coefficients  for  calculating  the  longitudinal  component 

VX 23.  VX24  of  velocity  of  the  inner  track  (VX2). 

VX2P  =  Array  for  storing  VX2. 

VXD  =  Derivative  with  respect  to  time  of  VX. 


VXP  =  Array  for  storing  VX. 

VXPRV  =  Previous  value  of  VX. 

VY  =  Lateral  component  of  velocity  of  the  vehicle,  in/sec. 
VYD  =  Derivative  with  respect  to  time  of  VY. 

VYP  =  Array  for  storing  VY. 

VYPRV  «  Previous  value  of  VY. 

W  =  Yaw  angle,  rad. 

WD  =  Yaw  rate,  rad/sec. 

WDP  =  Array  for  storing  WD. 

WDPRV  =  Previous  value  of  WD. 

WP  =  Array  for  storing  W. 

Z1  =  Slip  radius  of  the  outer  track. 

ZIP  =  Previous  value  of  Z1 . 

Z2  =  Slip  radius  of  the  inner  track. 

Z2P  =  Previous  value  of  Z2. 

B.4  SUBROUTINE  AGIL 

DIB  =  Shearing  deformation  of  soil  under  the  outer  track 

D2B  =  Shearing  deformation  of  soil  under  the  inner  track 

E  =  Steering  ratio. 

FCX  =  Longitudinal  component  of  inertial  force. 

FCY  =  Lateral  component  of  inertial  force. 

P  =  Offset. 

PT1  =  Power  required  by  the  sprocket  of  the  outer  track. 

PT2  =  Power  required  by  the  sprocket  of  the  inner  track. 

PTE  =  Power  required  from  the  engine. 

PTS  =  Total  power  required  by  the  sprockets. 

Q1  =  Lateral  component  of  shear  stress  along  the  outer  track 

Q2  =  Lateral  component  of  shear  stress  along  the  inner  track 

RIB  =  Normal  stress  under  the  outer  track. 

R2B  =  Normal  stress  under  the  inner  track. 

RIOL  =  Instantaneous  radius  of  curvature. 

ROL  =  Ordinate  of  the  instantaneous  center  of  rotation  of  the 
vehicle. 

SF  =  Coefficient  of  rolling  resistance. 
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T1  =  Longitudinal  component  of  shear  stress  along  the  outer 
track. 

T2  =  Longitudinal  component  of  shear  stress  along  the  inner 
track. 

TAU1  =  Maximum  shear  strength  of  soil  under  the  outer  track. 

TAU2  =  Maximum  shear  strength  of  soil  under  the  inner  track. 

TQ1  =  Shear  stress  under  the  outer  track. 

TQ2  =  Shear  stress  under  the  inner  track. 

V  =  Velocity  of  the  vehicle. 

VSX1  =  Longitudinal  component  of  slip  velocity  of  the  outer 
track . 

VSX2  =  Longitudinal  component  of  slip  velocity  of  the  inner 
track . 

VX  =  Longitudinal  component  of  velocity  of  the  vehicle. 

VX1  =  Longitudinal  component  of  velocity  of  the  outer  track. 

VX2  =  Longitudinal  component  of  velocity  of  the  inner  track. 

VXD  =  Time  derivative  of  VX. 

VY  =  Lateral  component  of  velocity  of  the  vehicle. 

VYD  =  Time  derivative  of  VY. 

W  =  Yaw  angle. 

WD  =  Yaw  rate. 

Z1  =  Slip  radius  of  the  outer  track. 

Z2  =  Slip  radius  of  the  inner  track. 

B.5  SUBROUTINE  MIVCI 

AMI  =  Mobility  index. 

BF  =  Bogie  factor. 

CF  =  Clearance  factor. 

CPF  =  Contact  pressure  factor. 

D  =  Track  width,  in. 

DRAD  =  Factor  for  converting  from  degrees  to  radians. 

EF  =  Engine  factor. 

GF  =  Grouser  factor. 

GH  =  Height  of  the  grouser,  in. 

HP  =  Specified  engine  horsepower,  hp. 

HPT  =  Horsepower  per  ton. 


L  -  Length  of  track  in  contact  with  the  ground,  in. 

NB  -  Total  number  of  bogies  on  track  in  contact  with  ground 
THETAA  =  Approach  angle  of  the  track  envelope,  deg. 

TRKF  «  Track  factor. 

TRNF  -  Transmission  factor. 

TRNT  =  Transmission  type  (0  for  manual,  1  for  automatic). 

TSL  »  Track  shoe  length,  in. 

VCI1  -  Vehicle  cone  index  for  one  pass. 

WF  =  Weight  factor. 

WT  =  Weight  of  the  vehicle,  lb. 

B.6  SUBROUTINE  CONE  (INPUT  AND  OUTPUT  ONLY) 

C  =  Static  cohesive  component  of  shear  strength,  psi. 

Cl  =  WES  cone  index. 

CL  -  Cone  length,  in. 

DI  =  Cone  diameter,  in. 

G  -  Shear  modulus,  psi/in. 

GAMA  =  Density,  psi/in. 

PHI  =  Angle  of  internal  friction,  deg. 

Z  «  Depth  of  penetration,  in. 


APPENDIX  C:  DESCRIPTION  OF  THE  FIELD  DIRECT  SHEAR  DEVICE 
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C.1  BACKGROUND 

The  terrain-vehicle  interaction  model  described  in  Chapter  2  requires  six 
soil  parameters  as  input  that  have  to  be  determined  experimentally.  The 
parameters  are  (Figures  2. 1-2. 3): 
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y 
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G  Initial  shear  stiffness  coefficient  (assumed  to  be  independent  of 
rate  of  deformation) 

A  Material  constant  describing  the  maximum  shearing  strength  of  the 
material  at  very  high  normal  loading  (Equation  2.1) 


M  Material  constant  related  to  the  parameter  A  and  to  the  static  soil 
cohesion  C  as  M  =  A  -  C 

N  Material  parameter  which  appears  in  Equation  2.1 

Cd  Increase  in  soil  cohesion  due  to  dynamic  loading  (maximum  value 

achieved  for  loading  rates  of  interest) 

A  Material  constant  describing  the  effects  of  rate  of  deformation  on 
the  cohesive  strength 

<p  Angle  of  internal  friction  if  Equation  2.2  is  used 
The  soil  parameters  G  ,  A  ,  M  ,  N  ,  and  $  can  be  determined  from  various 
existing  laboratory  test  devices,  such  as  the  triaxial  shear  device  or  direct 
shear  device.  The  triaxial  shear  and  direct  shear  devices,  however,  may  not 


yield  the  same  values  of  G  ,  A  ,  M  ,  N  ,  and  <j>  for  identical  specimens 


because  of  differences  in  test  boundary  conditions.  The  stress  boundary 
conditions  associated  with  the  direct  shear  test  more  closely  approximate  the 
stress  conditions  experienced  by  the  soil  during  steering  of  track-laying 


vehicles.  It  is,  therefore,  more  appropriate  to  determine  these  parameters 


from  direct  shear  tests.  The  parameters  Cd  and  A  can  only  be  determined 
from  special  static  and  dynamic  triaxial  shear  tests  since  dynamic  direct 


shear  devices  are  not  presently  available.  Therefore,  to  adequately  determine 


the  five  soil  parameters,  two  separate  test  series  may  be  required: 
1 .  Direct  shear  tests  to  define  G  ,  A  ,  M  ,  N  ,  and  p  . 


2.  Static  and  dynamic  triaxial  shear  tests  to  define  Cd  and  A  . 

It  should  be  noted  that  in  determining  Cd  and  A  from  triaxial  tests  rather 
than  direct  shear  tests,  it  is  assumed  that  these  parameters  are  not  sensitive 


to  test  boundary  conditions.  The  validity  of  this  assumption  should,  however, 


be  evaluated. 
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The  most  important  consideration  in  conducting  laboratory  soil  tests  is 
that  the  undisturbed  specimens  be  representati ve  of  the  materials  over  which 
the  vehicle  must  travel.  This  fact  implies  that  the  upper  several  inches  of 
surface  material  must  be  sampled,  trimmed  to  necessary  specimen  size,  and 
tested  in  the  laboratory.  Water  content,  soil  structure,  density,  and 
vegetation  root  systems,  all  of  which  affect  material  response,  must  be 
preserved.  With  this  in  mind,  a  field-operated  direct  shear  device  capable  of 
testing  a  variety  of  in  situ  surface  soils  for  normal  loads  of  interest  was 
designed  and  fabricated.  The  description  of  the  device  and  the  procedure  by 
which  the  soil  parameters  can  be  determined  are  documented  in  this  appendix. 

C . 2  DIRECT  SHEAR  DEVICE 

C.2.1  Design  Consideration 

Previously  proposed  field  devices  were  considered  but  rejected  because  of 
one  or  more  of  the  following  reasons:  some  of  the  soil  parameters  could  not 
be  measured  and  hence  required  additional  tests;  the  necessary  support 
equipment  was  too  massive  to  be  easily  field  transportable;  or  specimen 
disturbances  were  encountered  before  testing.  The  idea  of  creating  a  new  type 
of  test  was  also  rejected  because  any  new  device  would  contain  inherent 
boundary  problems,  all  of  which  would  have  to  be  evaluated  with  time  and 
usage.  The  direct  shear  device,  on  the  other  hand,  has  been  used  extensively, 
and  the  direct  shear  test  is  a  simple  test  to  run.  Furthermore,  the  five 
basic  soil  parameters  (G  ,  A  ,  M  ,  N  ,  and  ((>)  could  be  measured  directly 
from  this  test.  Figure  C.1  shows  a  sketch  of  the  field  device  that  was 
fabricated  as  a  result  of  this  project.  Figure  C.2  contains  photographs  taken 
of  the  device  during  the  conduct  of  actual  field  tests. 

C.2. 2  Specimen  Container 

Specimen  configuration  was  the  first  consideration  made  in  the  design  of 
the  device.  It  was  assumed  that  in  many  cases  the  in  situ  soil  could  not  be 
sampled  without  disturbance;  therefore,  the  specimen  container  would  have  to 
be  placed  around  the  soil.  A  round  ring  similar  to  a  coring  device  would 
afford  the  least  chance  of  soil  disturbance.  However,  the  stress  distribution 
along  a  horizontal  plane  of  a  circular  specimen  is  not  uniform.  To  reduce  the 
nonuniformity,  a  square-shaped  specimen  container  was  selected. 
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A  4-  by  4- in  box  was  selected  in  order  to  keep  the  shear  and  normal  loads 

within  limits  of  interest  to  the  analysis  of  track-laying  vehicles  and  at  the 

same  time  retain  a  reasonably  large  specimen  size.  The  use  of  deadweights  is 

the  simplest  way  to  produce  normal  load,  but  use  of  more  than  200  lb  in 

weights  is  awkward  for  field  testing.  Therefore,  with  the  weight  requirement 

below  the  200-lb  limit,  normal  stress  of  up  to  12  psi  can  be  produced  on  a 

2 

4-  by  4-in,  or  16-in  specimen.  However,  the  largest  particle  or  grain  size 
permissible  with  a  4-  by  4-in  specimen  is  probably  1/2  in,  which  is  a 
reasonable  limit  for  most  terrains  of  interest. 

The  overall  specimen  height  was  controlled  by  the  depth  of  the  desired 
shear  plane  as  directed  by  grouser  depth  ranging  from  approximately  3/4  to 
1-1/2  in.  The  compressibility  of  soil  could  significantly  alter  this  depth, 
but  for  estimation  purposes  the  depth  was  assumed  to  be  no  greater  than 
2  in.  Therefore,  the  height  of  the  upper  box  portion  was  set  at  2  in, 
permitting  testing  of  depths  from  approximately  1/4  to  2  in.  This  height  can 
be  altered  should  particular  site  conditions  dictate.  The  lower  box  portion 
was  set  at  a  1-1/4-in  height,  including  the  cutting  edge.  A  1/8-in  wall 
thickness  was  used  for  both  boxes. 

Figure  C.3  presents  a  series  of  sketches  of  the  specimen  container  show¬ 


ing  the  various  stages  of  placement.  To  minimize  specimen  disturbance,  it  was 
decided  to  use  the  specimen  cutting  box  as  the  device  container  rather  than 
remove  the  cutting  and  place  a  container  over  the  specimen.  The  box  consists 
of  three  parts:  (a)  a  lower  portion  with  knife-sharp  edges  to  aid  in  cutting 


the  soil,  (b)  an  upper  portion,  and  (c)  an  outer  holder  to  keep  the  lower  and 
upper  portions  in  alignment.  The  box  is  alternately  pushed  and  trimmed  into 
the  soil  to  the  desired  depth.  Once  in  place,  the  outer  holder  can  be 
carefully  removed,  leaving  the  two  boxes  on  the  specimen  with  the  joint 
between  the  lower  and  upper  box  portion  forming  the  shear  test  plane. 

C.2.3  Base 

A  relatively  narrow  1 -in-thick  aluminum  plate  was  used  to  construct  the 
base  with  a  square  hole  at  one  end  to  fit  around  the  4-  by  4-in  lower  specimen 
containers  (Figure  C.1).  The  shear  loader  was  attached  to  the  other  end  of 
the  plate.  A  second  1-in-thick  aluminum  yoke  was  constructed  to  fit  over  the 
upper  specimen  container.  Set  screws  through  the  yoke  serve  to  raise  the  yoke 
off  the  baseplate,  thus  minimizing  friction  between  the  surfaces.  The  shear 
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loader  attached  by  cable  to  the  yoke  pulls  the  upper  specimen  while  the  base 
reacts  against  the  lower  specimen  container.  Guide  rails  along  the  edge  of 
the  base  ensure  that  no  torsional  shear  deformation  or  twisting  is  applied  to 
the  specimen. 

C.2.JJ  Shear  Loader 

An  electric  12-volt  boat  winch  was  incorporated  into  the  base  as  the 
shear  loader.  This  approach  is  the  simplest  for  providing  a  shear  loader. 

(If  necessary,  the  winch  can  be  replaced  with  a  more  sophisticated  loader 
custom  built  for  this  device.)  Currently,  the  winch  is  capable  of  pulling 
loads  up  to  2000  lb.  Static  loading  rates  can  be  applied  by  manually  turning 
the  winch  via  a  socket-rachet  arrangement.  Fast  loading  rates  (approximately 
300-600  msec  time  to  peak  load  that  is  equivalent  to  a  strain  rate  of  0.5  to 
1.5  per  second,  which  is  compatible  with  the  strain  rate  under  the  track)  can 
be  applied  using  the  electric  feature  of  the  winch.  The  power  is  supplied  by 
a  12-volt  car  battery,  which  is  also  used  as  the  instrumentation  power  supply. 

C.2.5  Instrumentation 

A  2-in  travel  film  potentiometer  is  attached  to  the  base  and  records 
relative  movement  between  the  upper  specimen  holder  and  the  base.  A  strain 
gage  load  cell  attaching  the  winch  cable  to  the  specimen  yoke  is  used  to 
measure  shear  load.  A  compact,  two-channel  DC  instrumentation  amplifier  is 
used  for  signal  conditioning.  Output  is  recorded  in  the  form  of  a  shear  load 
versus  deflection  plot  on  a  commercially  available  DC-operated  X-Y  plotter. 

As  previously  mentioned,  a  simple  car  battery  is  the  main  power  supply.  All 
initial  testing  was  done  by  recording  the  data  on  a  time  base  light  beam  strip 
chart.  This  recording  procedure  was  later  dropped  since  the  loading  times 
remained  fairly  constant  on  the  soils  tested.  A  time  base  can  be  added  at  a 
later  date  through  the  use  of  a  frequency  oscillator  and  an  X-Y-Z  recorder. 

C.2.6  Normal  Load 

A  series  of  steel  weights,  the  largest  weighing  57  lb  and  the  smallest 
weighing  8-1/2  lb,  was  fabricated  for  use  with  the  device.  Guide  holes  and 
studs  permit  stacking  and  centering  of  the  weights  on  the  specimen  top 
surface.  Although  a  variety  of  load  combinations  are  possible,  most  tests 
have  been  conducted  using  weights  totaling  approximately  8.6,  36.6,  65.7,  and 
122.7  lb  (i.e.,  normal  stress  levels  of  0.51*,  2.29,  4.11,  and  7.67  psi). 


C.3  MEASUREMENT  OF  SOIL  PARAMETERS 

A  series  of  two  or  more  tests  is  required  at  a  site  to  define  the 
necessary  soil  parameters.  A  typical  testing  program  may  call  for  the  conduct 
of  four  fast  and  four  slow  tests  at  normal  stresses  of  0.54,  2.29,  4.11,  and 
7.67  psi.  For  each  test,  an  X-Y  data  record  of  shear  load  versus  deflection 
is  obtained.  In  addition,  measurements  of  soil  density  and  water  content  are 
made  on  each  test  specimen  (generally  on  the  posttest  specimen  contained  in 
the  upper  and  lower  specimen  holders). 

For  each  test,  a  plot  of  shear  stress  versus  deflection  is  obtained.  The 
initial  slope  of  the  plot  defines  G  ,  the  peak  stress  defines  the  maximum 
shear  stress,  and  the  deflection  at  peak  stress  divided  by  time  to  peak  stress 
defines  the  deflection  rate.  A  table  listing  of  each  test  is  used  to  summarize 
the  data  and  contains  the  specimen  number,  wet  density,  water  content,  dry 
density,  normal  load/stress,  maximum  shear  load/stress,  initial  G  , 
deflection  at  peak  stress,  and  deflection  rate.  Figure  C.4  presents  the  test 
results  obtained  from  the  series  of  field  tests  conducted  at  a  given  site. 

The  analysis  plots  are  shown  graphically  in  Figure  C.5.  A  summary  plot 
of  3hear  stress  versus  shear  deformation  is  made  to  obtain  either  static  or 
dynamic  failure  envelopes.  From  the  static  failure  envelope,  the  values  of 
A  ,  M  ,  N  ,  and  <p  are  obtained.  The  value  of  Cd  and  A  are  obtained 
from  the  dynamic  failure  envelope  as  shown  in  the  plot  of  cohesion  versus  rate 
of  deformation.  The  value  of  G  is  the  initial  slope  of  shear  stress-shear 
deformation  curve  (Figure  C.4). 
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Figure  C.3.  Cross  sections  through  the  specimen  container 
showing  various  stages  of  placement. 
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Figure  C.5.  Graphical  presentation  of  the  selection  of 
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APPENDIX  D:  NOTATION 


Forward  acceleration  of  the  vehicle 
Lateral  acceleration  of  the  vehicle 
Acceleration  of  the  vehicle  along  the  4>  axis 
Acceleration  of  the  vehicle  along  the  y  axis 
Material  constant  in  failure  envelope 


Track  tread  (distance  between  the  centerlines  of  the 
tracks) 

CL2/W 

CdL2/W 


Static  cohesive  component  of  shear  strength 

Added  cohesive  strength  due  to  dynamic  loading 

Abscissa  of  the  center  of  gravity  of  the  vehicle 

Slip  radius  of  the  outer  track 

Slip  radius  of  the  inner  track 

Center  of  gravity  of  the  vehicle 

WES  cone  index 

Center  of  rotation  of  the  vehicle 


Track  width 
Cone  diameter 

FCX/W 

fcy/w 

Inertial  force 

Longitudinal  component  of  inertial  force 
Transverse  component  of  inertial  force 
Coefficient  of  rolling  resistance 
Acceleration  due  to  gravity 
Initial  shear  stiffness  coefficient 


Height  of  center  of  gravity 


Mass  moment  of  inertia  of  the  vehicle  about  an  axis 
passing  through  its  center  of  gravity  and  parallel  to 
the  2  axis 

Center  of  slip  rotation  of  the  outer  track 

Center  of  slip  rotation  of  the  inner  track 

Instantaneous  center  of  rotation  of  the  vehicle 

Distance  between  adjacent  wheels 

Length  of  the  track  in  contact  with  the  ground 

Cone  length 

ML2/W 

Material  constant  in  failure  envelope 

Mobility  index 

NW/L2 

Vertical  component  of  T 

Vertical  component  of  T2 

Material  constant  in  failure  envelope 

Contribution  due  to  the  outer  track  tension 

Contribution  due  to  the  inner  track  tension 

P/L 

Offset  (distance  from  center  of  gravity  to  pivot  point 
of  vehicle) 

Total  power  =  PT1  +  PT2 

Power  required  by  the  sprocket  of  the  outer  track 

Power  required  by  the  sprocket  of  the  inner  track 

Differential  power  =  PT1  -  PT2 

Power  required  from  the  engine 

dL2Q1 (x)/W 

dL2Q2(x)/W 

Transverse  component  of  shear  stress  along  the  outer 
track 

Transverse  component  of  shear  stress  along  the  inner 


track 

ar  2d 


R,  (X) 
R2(X) 


(x) 
t2(x ) 


T,  (X) 


T2(X) 


,v  ,  v 
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Ordinate  of  the  instantaneous  center  of  rotation  of  the 
vehicle 

Radius  of  trajectory  of  center  of  gravity  of  vehicle 

Rolling  resistance 

Normal  stress  under  the  outer  track 

Normal  stress  under  the  inner  track 

Instantaneous  radius  of  curvature 


dLrT  1  (x  )/W 
dL2T2(x)/W 

Track  tension  in  the  inner  track 
Track  tension  in  the  outer  track 

Longitudinal  component  of  shear  stress  along  the  outer 
track 

Longitudinal  component  of  shear  stress  along  the  inner 
track 

Velocity  of  the  vehicle 

Instantaneous  velocity  of  an  arbitrary  point  of  the  hull 

and  its  components  along  X  and  Y  coordinates 

Total  slip  velocity  of  the  outer  track 

Total  slip  velocity  of  the  inner  track 

Longitudinal  component  of  slip  velocity  of  the  outer 

track 

Longitudinal  component  of  slip  velocity  of  the  inner 
track 

Transverse  component  of  slip  velocity  of  the  outer  track 
Transverse  component  of  slip  velocity  of  the  inner  track 
Longitudinal  component  of  velocity  of  the  vehicle 
Longitudinal  component  of  velocity  of  the  outer  track 
Longitudinal  component  of  velocity  of  the  inner  track 
Transverse  component  of  velocity  of  the  vehicle 
Component  of  velocity  of  the  vehicle  along  the  $  axis 
Component  of  velocity  of  the  vehicle  along  the  T  axis 
Vehicle  cone  index  for  one  pass 
Weight  of  the  vehicle 


Wjj  Component  of  weight  of  the  vehicle  normal  to  the  terrai 

WT  Component  of  weight  of  the  vehicle  parallel  to  the 

terrain 
x  X/L 


X.Y.Z 

y 


2 

z 

a 

a 
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Y 


A 


e 


n 


e 


A 

x 
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Local  coordinate  system 

Y/L 

Z/L 

Depth  of  penetration  of  the  cone 
Side-slip  angle 

Half  of  the  apex  angle  of  the  cone 
l/L 

Density  of  the  material 

Angle  of  slip  direction  of  the  outer  track 
Angle  of  slip  direction  of  the  inner  track 
Shearing  deformation 

Initial  displacement  of  the  outer  track 

Initial  displacement  of  the  inner  track 

Shearing  deformation  of  soil  under  the  outer  track 

Time  rate  of  shearing  deformation 

Shearing  deformation  of  soil  under  the  inner  track 

Time  rate  of  shearing  deformation 

A,/L 

A/L 

a2/l 

a2/l 

Steering  ratio 

Angle  of  sloping  terrain 

Directional  angle 

Approach  angle  of  the  track  envelope 
Departure  angle  of  the  track  envelope 
Material  constant  related  to  rate  effect 
AL 

GL3/W 


v 


Angle  related  to  initial  position  of  the  vehicle 


0 


Normal  stress 
x  Shear  stress 

Maximum  shear  strength 
0  Angle  of  internal  friction 

Coordinate  system  fixed  on  level  ground 

u>  Yaw  angle 
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